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THE  WATERS  OF  HOT  SPRINGS  NATIONAL  PARK,  ARKANSAS— 
THEIR  ORIGIN,  NATURE,  AND  MANAGEMENT 

By  M.  S.  Bedinger,  F.  J.  Pearson,  Jr.,  J.  E.  Reed, 
R.  T.  Sniegocki,  and  C.  G.  Stone 

ABSTRACT 

The  47  hot  springs  of  Hot  Springs  National   Park,  Ark.,   issue  from 
the  plunging  crestline  of  a   large  overturned  anticline,  along  the 
southern  margin  of  the  Ouachita  anti cl inorium,   in   the  Zigzag  Mountains. 
The  rocks   in   the  vicinity  of  the  hot  springs   range  in  age  from  Ordovician 
to  Mississippian.     The  rocks— cherts ,  novaculites,  sandstones,  and 
shales— are  well    indurated,   folded,   faulted,   and  jointed.     The  springs 
emerge   from  the  Hot  Springs  Sandstone  near  the  anticlinal    axis,  between 
the  traces  of  two  thrust  faults   that  are  parallel   to  the  axis   of  the 
anti  cline. 

The  combined  flow  of  the  hot  springs   ranges   from  750,000  to 
950,000  gallons  per  day   (3.28X10"^   to  4.16X10"^   cubic  meters   per  second). 
The  flow  of  the  springs   fluctuates  seasonally,  corresponding  to  the 
annual   climatic  cycle.     Flow  is  highest  in  the  winter  and  spring  and  is 
lowest  in   the  summer  and  fall.     The  temperature  of  the  combined  hot- 
springs  water  is  about  62°  Celsius. 

The  radioactivity  and  chemical  concentration  of  the  hot  springs  are 
similar  to  that  of  cold-water  springs  and  wells   in  the  area.     The 
dissolved-solids  concentrations  of  the  waters  in  the  area  generally 
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range   from   175   to  200  milligrams  per  liter.      The  main   differences    in  the 
quality  of  the  hot  water,   compared  with  nearby  cold  ground  waters,  are 
the  high   temperature  and  the  higher  silica  concentration  of  the  hot  springs. 
Cold  waters   in   the  area  generally  range  from  15.0°   to  26.8°  Celsius. 
The  silica  concentrations  of  cold  ground  waters  range  from  2.6  to  13.0 
milligrams   per   liter,  whereas   silica  concentration  of  the  hot  springs   is 
about  40  milligrams   per   liter.     The  high  silica  concentration  of  the  hot 
springs   is  due  to  the  increased  solubility  of  silica  in  hot  water.     The 
silica  concentration  of  the  hot  springs   indicates  that  the  maximum  tempera- 
ture reached  by  the  hot-springs  water  is  no  more  than  a  few  degrees  higher 
than   the   temperature  at  which   the  springs  emerge. 

The  tritium  and   ^"^C  analyses  of  the  water  indicate  that  the  water 
is   a  mixture  of  a  very  small   amount  of  water  less   than  20  years  old 
and  a  preponderance  of  water  about  4,400  years  old.     The  deuterium  and 
^^0  concentrations  of  the  cold  ground  waters  are  not  significantly  different 
from  that  of  the   hot  spring  water. 

The  presence  of  radium  and  radon  in  the  hot-springs  waters   has 
been  established  by  analyses.     Recent  (1973)   analysis  showed  the  radium 
concentration  to  be  2.1   picocuries  per  liter.     Analyses  of  the  radon  gas, 
a  radioactive  decay  product  of  radium,  made  in  1953,   ranged  from  0.14  to 
30.5  nanocuries   per  liter. 

The  geochemical    data,   flow  measurements,  and  geologic  structure 
of   the   region   support   the  concept   that  virtually   all    the  hot-springs 
water  is    local,   meteoric  water.      Recharge   to  the  hot-springs   artesian- 
flow  system  is   by   infiltration  of  rainfall    in   the  outcrop  areas  of  the 


Bigfork  Chert  and   the  Arkansas  Novaculite.     The  water  moves   slowly   to 
depth  where  it  is   heated  by  contact  with   rocks  of  high   temperature. 
Highly  permeable  zones,   related  to  jointing  or  faulting,  collect  the 
heated  water  in   the  aquifer  and  provide  avenues   for  the  water  to  travel 
rapidly  to  the  surface. 

INTRODUCTION 
0 b je^ tij/ e s_and  Scope  of  the  Study 

The   thermal   springs  of  the  Hot  Springs  National    Park,  at  Hot 
Springs,  Ark.    (fia.    1),   have  been  a  natural    resource  of  international 
reno'./n   for  many  years.      The  springs  were   known   to  President  Thomas 
Jefferson,  who   initiated   the   first  scientific  study    in   1804.      This   study, 
by  William  Dunbar  and  George  Hunter,  marked  the  beginning  of  an  era  of 
scientific  curiosity  as    to   the  origin  and  heat  source  of  the  springs. 
With   the   advent  of  modern-day  environmental    problems  —  effects   of  un- 
controlled urbanization   and  overuse  of   the  natural    resources — scientific 
knowledge  of  the  springs   becomes   necessary   for  management  of  the  environ- 
ment and  the   resources    for  beneficial    use  by   future  generations.     It  is 
in   this   context  that  this  study  was  made  by  the  U.S.   Geological   Survey, 
on  behalf  of  the  National    Park  Service,    to  investigate  the  origin  of 
the  thermal  water  resources   and  the  geohydrologic  framework  of  the  Hot 
Springs   National    Park. 

The  purpose  of  this  study   is    to  (1)   describe  the  hydrogeologic 
f'-ampwork  -ind   functioninq  of  the  spring  system;   (2)   consider,   in  accordance 
with   the  most  acceptable  hypothesis  of  spring  origin,   the  possible  effects 
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Figure  1.-  Location  of  study  area. 


of  changed  land  use  and  urban  development  on   the  hot  springs;  and 
(3)    recommend  an  investigational   and  monitoring  program  that  will    insure 
adequate   records   on  which  to  base  management  and  operational    decisions 
for  best  use  of  the  spring  waters   and  preservation  of  esthetic  values. 

The  process  by  v/hich  water  enters   the  spring  system,  becomes 
heated,   and  discharges  at   the  hot  springs  has  been   identified  using 
updated  geologic  mapping,  modern   knowledge  of  ground-water  hydrology, 
application  of  geochemical    and  isotopic  studies,   and  data  on  volume, 
temperature,  and  chemical   quality  of  the  spring  discharge.     Old  records 
have  been  examined  to  determine  their  applicability  to  the  current 
problem,   and  as   a  basis   for  planning  data  collection   for  the  present 
study  and  for  planning  continuing  data  collection  for  future  management 
of  the  springs. 

Public   interest   in   the  hot  springs   has  been  focused  primarily  on 
the  therapeutic  value  of  the  waters,   and,   in  serving  such  interest,  this 
also  has  been  the   focus  of  Federal   management  since  the  area  was  estab- 
lished as   the  Hot  Springs  Reservation,   in   1832.     This  emphasis   did  not 
change  when  direct  Federal   supervision  was   implemented,   in  1877,  and 
when  the  area  was  designated  as   a  National   Park,   in   1921.     The  purpose 
of  the  park   today   is   to  preserve  and  protect  the  hot  springs   for  present 
and  future  generations. 

Long-range  planning  for  park  values   and  uses   takes   into  consideration, 
however,   the  prospect  of  a  shifting  of  emphasis  from  therapeutic  values 
of  the  spring  waters   to  the  scientific,  esthetic,  and  recreational 
values  of  the  park  as   a  whole.     The  probable   life  of  the  springs  and  their 


thermal  characteristics  play  an  important  role  in  attracting  visitors 
to  the  area.  For  this  reason,  this  report  and  the  investigational 
program  recommended  herein  emphasize  the  origin  and  functioning  of  the 
spring  system  and  the  possible  effects  of  man-imposed  stresses  on  the 
future  of  the  springs  and  the  park. 

Metric  Units 

For  those  readers  interested  in  using  the  metric  system,  metric 

equivalents  of  English  units  are  given  in  parentheses.  The  English 

units  in  this  report  may  be  converted  to  metric  units  by  making  the 

following  conversions: 

To  convert  from       Multiply  by  To  obtain 

(English  unit)     (conversion  factor)      (metric  unit) 

Foot  (ft)  3.048X10"^  meter  (m) 

Square  mile  (mi-)  2.59X10  square  meter  (m^) 

Gallons  per  day  (gpd)  4.38X10"  cubic  meters  per  second  (m^/s) 

Feet  per  day   (ft/day)  3.53X10"  meters   per  second  (m/s) 

The  conversion   from  temperature   in   degrees   Fahrenheit  (°F)   to 
temperature   in  degrees   Celsius   (°C)    is  expressed  by   :     °C=(5/9) (°F-32). 
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History  of  Spring  Development  and  Scientific  Study 

The  history  of  the  hot  springs  area  has  been  documented  in  numerous 
publications,  many  of  which  present  detailed  accounts  of  some  aspect  of 
the  springs'  environment  and  the  cultural  development  of  the  area.  For 
the  purposes  of  this  report,  therefore,  and  to  minimize  duplication, 
only  those  historic  events  and  developmental  practices  that  relate  to 
the  technical  management  of  the  springs  are  cited. 

Early  descriptions  of  the  hot  springs  give  different  accounts  of 
as  many  as  72  spring  openings,  in  a  belt  about  one- fourth  mile  long 


and  a  few  hundred  feet  wide,  along  the  southwest  slope  of  Hot  Springs 
Mountain.     Excavation  and  covering  of  springs,   to  increase  and  concen- 
trate flows  and  to  protect  the  springs   from  contamination,  have  so 
altered  the  natural   spring  environment  that  it  bears   no  resemblance 
to  the  original   condition.     Among  the  early  investigators,  Owen   (1860) 
reported  42  springs;  Glasgow  (1860),   54  springs;  Haywood  (1902),  46 
springs;  and  Hamilton   (1932),  48  springs.      In  his   detailed  history  of 
Hot  Springs,  Scully  (1966,  p.    139)    reported  47  active  springs,   including 
2  exhibition  springs. 

Prior  to  1877  some  of  the  springs  were  walled  in  and  covered  by 
masonry  arches   to  protect  them  from  contamination   (Scully,   1966,  p. 
118).     By  the  1890's,  most  of  the  springs  were  covered  and  a  complicated 
piping  system  had  evolved  for  supplying  the  bathhouses  with  hot  water. 
In   1901    the  springs  were  uncovered  to  give  access   for  sampling,  and 
chemical   analyses  were  made  by  J.   K.   Haywood  (1902).     The  springs  were 
opened  again  in  1931   for  cleaning;  some  of  them  were  deepened,  and  the 
present-day   (1974)   collection  system  was  constructed.     The  collecting 
system  diverts   the  flow  of  44  springs   to  a  central    reservoir,   from  which 
the  water  is    redistributed  to   individual   bathhouses.      A  map  of  the  col- 
lection system  is  shown   in  figure  2.     Since  1948  all    the  water  delivered 
to  the  bathhouses  has  been  metered.     Excess  water  overflows  into  Hot 
Springs   Creek  when  storage  reservoirs   are  full. 

Through  the  years   (1860  to  the  present),  at  least  20  scientific 
investigations,  directly  or  indirectly  involving  the  hot  springs,  have 
been  made.     Although  each  study  generally  had  a  separate  and  specific 
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objective,  many  of  the  investigators  became  sufficiently  interested  in 
the  hot  springs  to  try  to  explain  the  origin  of  the  water  and  the  source 
of  the  heat. 

The  chemical  quality  of  water  from  the  hot  springs  in  Arkansas  has 
been  of  great  interest  to  man,  probably  since  the  hot  springs  were  first 
discovered.  One  of  the  earliest  scientific  approaches  to  determine  the 
concentration  of  the  hot  springs  waters  is  found  in  J.  C.  Branner's  (1892) 
Annual  Report  for  1891,  in  which  analyses  of  water  samples  collected  in 
1890  dre   tabulated  in  grains  per  gallon.   At  random  times  since  1890, 
analyses  have  been  made  for  i nvestiaati ons.  The  purpose  of  many  of 
these  investigations  has  been  to  support  some  theraputic  claim  for  the 
water,  or  to  determine  whether  the  chemical  concentration  of  the  water 
has  changed. 

Most  earlier  investigators  concluded  that  the  waters  discharged 
from  the  hot  springs  are  of  meteoric  origin,  having  fallen  as  precipita- 
tion and  recharged  to  the  Bigfork  Chert  in  the  anticlinal  valley  lying 
just  northwest,  north,  and  northeast  of  Hot  Springs.  Some  investigators 
have  attributed  some  of  the  recharae  to  the  outcrop  of  Arkansas  Novaculite 
to  the  east  of  the  hot  springs.  Another  theory  that  is  regarded  by  some 
to  have  a  degree  of  scientific  validity  is  that  the  water  may  be  of 
juvenile  origin,  that  is,  derived  from  the  interior  of  the  earth  and  not 
having  previously  existed  as  atmospheric  water. 

Bryan  (1922,  p.  426)  posed  the  question  as  to  the  meteoric,  juvenile, 
or  mixed  origin  of  the  waters  discharged  from  the  hot  springs.  He  indicated 
(p.  447-448)  that  the  juvenile  theory  is  perhaps  more  satisfactory,  although 


it  rests   on  an   insecure  foundation  in  postulating  (1)   a  special    igneous 
mass   that  is  discharging  water  owing  to  cooling  and  recrystal lization, 
and  (2)   a  special   fault  fissure   through  which  the  water  rises   to  the 
land  surface.     Bryan   (p.   444)   analyzed  the  merits  of  both  the  juvenile 
and  the  meteoric  theories,  but  conceded  that  ".    .    .a  definite  con- 
clusion as   to  the  ultimate  origin  of  the  water  in  the  Hot  Springs   cannot 
now  be  reached."     He  pointed  out  (p.    443-444)   that  "If  the  water  is 
juvenile  there  is  presumably  a  constant  supply,  diminishing  ^^ery  gradually 
through  the  centuries   in  quantity  and  temperature   .    .    .   If  .    .    .   the  water 
has   a  meteoric  origin,   it  is  variable  in  quantity,  fluctuating  with  the 
seasons  or  with  groups  of  years  having  heavy  or  light  rainfall."     Thus, 
Bryan   recognized  the  critical   value  of  precise  measurements  of  temperature, 
discharge,  and  other  parameters  during  a  sufficient  period  of  time  to 
provide  adequate  data  on  which  to  base  conclusions   as   to  the  water's 
origin. 

Arndt  and  Stroud  (1953)  suggested  a  dual   origin  for  the  water. 
Meteoric  water,   they  believed,  entered  the  spring  system  through  the 
lower  member  of  the  Arkansas  Novaculite  on  Hot  Springs   and  North  Mountains. 
They  calculated  that  this  source  of  meteoric  water  could  supply  about 
one-sixth  of  the  flow  of  the  springs.     The  rest  of  the  water,  they  con- 
sidered,  could  be  juvenile  water  rising   from  depth. 

Proponents  of  the  theory  of  meteoric  origin  of  the  spring  waters 
include  Weed  (1902),  Purdue  (1910),   and  Purdue  and  Miser  (1923).     Purdue 
(1910,  p.   283)   described  the  geologic  conditions  supporting  this   view, 
and  identified  the  collecting  area  as  the  anticlinal   valley  between 
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Sugarloaf  and  West  Mountains.     Most  of  the  valley   is   underlain  by  the 
Bigfork  Chert,  a  much  fractured  formation  of  high  permeability  into 
which  infiltrates   a  large  part  of  water  from  precipitation.     According 
to  Purdue  (p.    284),   the  occurrence  of  this   formation   in  anticlinal 
valleys,  with  its  highly  inclined  beds,   affords  the  most  favorable  con- 
dition for  the  intake  of  water.     He  postulates   further  that  the  water 
passes   through   the  Bigfork  Chert,  beneath  the  North  Mountain  syacline, 
and  is   forced  upward  into  the  Hot  Sprinns   anticline  to  emerge  at  the 
hot  springs. 

This   suggested  movement  of  water  from  the  recharge  area  required 
geologic   conditions    that  account  for  passage  of  the  water  through  the 
Polk  Creek  Shale,  Missouri   Mountain  Shale,  and  the  Arkansas  Novaculite, 
to  discharge  as    it  does   from  the  Hot  Springs  Sandstone.     Such  conditions 
would  ordinarily   require  a  fault  or  faults,  with  associated  jointing 
and  fissuring,   that  would  provide  passage  through  these  formations. 
Several   authors  have  shown  such   a  fault  (Arndt  and  Stroud,   1953; 
Fellows,   1966).     Recent  mapping  by  Haley  and  Stone  (fig.    3)   confirms 
the  presence  of  a  complex  fault  system  in   the  area  but  indicates  no 
conclusive  evidence  of  a   large  fault  at  the  hot  springs.      In  addition 
to   the   faulting,   the   intensive  folding  and  overturning  of  formations   in 
the  vicinity  of  Hot  Springs  are  attended  by  intensive  jointing  and  fis- 
suring of  the  sandstone,   chert,  slate,   shale,   and  novaculite.     It  is 
concluded  from  the  present  study   that  the   faults  and  the  associated 
joints   and  fissures   provide  conduits   for  the  water. 
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Prof.   D.    D.   Owen,   in  his    report  of  1860  on  the  hot  springs,   said 
the  following:     "When  we  reflect  on   the  boundless  and  never-ceasing 
flow  of  thermal  waters   that  must  have  bathed  the  sides   of  Hot  Springs 
Ridge  for  countless  ages    .    .    .   and  however  inexplicable  such  wonderful 
phenomena  and  changes  may  at  first  appear,  yet,  when   the  chemical 
principles   become  properly   understood,   disclosed  by   the  enlightened 
and  accurate  chemical    analyses,   these  obscure  geological   transformations 
[and  the  origin  of  the  water  and  operation  of  the  springs]  can  be 
satisfactorily  and  clearly  explained."     It  would  seem  that  Owen's 
prophecy  may  finally  be   fulfilled. 

The  existence  of  juvenile  water  elsewhere  in   the  world  is  not  in 
question — water  vapor  has  been  observed  emanating  from  volcanoes   in 
dryland  areas  where  the  rocks   of  the  earth's   crust  contain   little  or  no 
meteoric  water.     However,   there  is  no  evidence  to  support  the  concept  of 
existing  large  quantities   of  juvenile  water  comparable  to  the  volumes  of 
water  known  to  be  of  meteoric  origin   (McGuiness,   1963). 

GEOLOGIC  SETTING 

The  rocks   cropping  out  in  the  vicinity  of  the  hot  springs  are  sedi- 
mentary  rocks,   although  intruded  igneous   rocks   are  exposed  in  the  region 
(Purdue  and  Miser,   1923).     The  sedimentary  rocks  are  relatively  old 
(Paleozoic)   and  consist  of  shale,  slate,   chert,  novaculite,  and  sandstone. 
The  names  of  the  geologic  formations,   their  geologic  ages,  and  positions 
in  the  geologic  column  are  given   in   table  1. 
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Though  no  significant  igneous   rocks  are  exposed  in  the  immediate 
vicinity  of  the  hot  springs,   their  nearby  occurrence  has  been  frequently 
cited  in  literature  as   possible  sources   for  the  heat  of  the  springs. 
The  igneous    rocks  were  intruded  into  the  sedimentary  rocks  during  the 
early  Late  Cretaceous   time  (about  90  million  years   ago).     The  larger 
igneous   intrusions   in   the  hot-springs   region  are  exposed  in   two  small 
areas,   about  6  miles  southeast  of  the  hot  springs.     Elsewhere  in  the 
region,  igneous   intrusions  occur  as   very  small   dikes   and  sills. 

The  sedimentary   rocks   in  the  vicinity  of  the  hot  springs  were 
originally  laid  down  on   the  bottom  of  the  sea  in  nearly  horizontal   beds. 
At  present  the  beds   are  generally  steeply  inclined,  because  of  tremendous 
mountain-building  forces   from  the  south  in   late  Paleozoic  time.     A 
geologic  map  (fig.    3)   shows   the  edges  of  the  inclined  strata  where  they 
intersect  the  land  surface.     When  the  formations   are  crossed  from  north- 
west to  southeast,   they  are  seen  in  cross  section   (fig.    3)   to  lie  in  a 
series  of  very  complexly  folded  anticlines   and  synclines,  with  some 
associated  thrust  faults. 

The  hot  springs  emerge  from  the  plunging  crestline  of  a  large  over- 
turned anticline  along  the  southern  margins  of  the  Ouachita  anticlinorium 
in   the  Zigzag  Mountains.     The  Zigzag  Mountains  basically  owe  their  presence 
to  the   resistant  exposures  of  the  Arkansas  Novaculite.     The  zigzag  pat- 
tern of  the  strata  is  mostly  due  to  tightly  compressed  folds  which  plunge 
southwestward  into  the  Mazarn  Basin.     The  Mazarn  Basin  is  a  structual   and 
topographic  basin   lying  south  of  the  Zigzag  Mountains.     The  Stanley  Shale, 
a  formation  much  less   resistant  to  erosion  than  the  Arkansas  Novaculite, 
crops   out  at  the  surface  of  the  Mazarn  Basin.     The  structual   setting  is 
illustrated  by  figure  6  in  Purdue  and  Miser  (1923). 
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The  formations,  composed  predominantly  of  shale,   include  the 
Womble,  Missouri  Mountain,  Polk  Creek,   and  Stanley  .     The  shales  have 
low  permeability,  but,   locally,   limestones   in   the  Womble  yield  water 
to  springs.     Shales  generally   impede  qround-water  movement,  except 
where  open  joints  and  fractures   are  present.     Wells   in  shales   generally 
yield  meager  quantities  of  water;   recharge  to  shales   is   also  small. 

The  Bigfork  Chert  typically   is   highly  permeable,  exhibiting  inter- 
granular  and  fracture  permeability.      The  Biqfork   is   composed  of  silt- 
sized,   generally  poorly  cemented  siliceous   particles    in   thin  beds   one- 
half  to  4  inches   (0.013-0.13  meter)    thick,  which  have  been  weathered, 
leaving  a  friable  material,   and  interbedded  with  lavers  of  dense  chert 
4  to   12   inches    (0.13-0.45  meter)   thick.     The  dense  chert  beds  were 
rendered  permeable  by  fracturing,  which  accompanied  the  intense  folding 
of  the  beds,  whereas   the  decalcified  silt-sized  material   has  significant 
intergranular  permeability  near  the  ground  surface. 

Wells   that  yield  the   largest  quantities   of  water  in   the  region  tap 
the  Bigfork  Chert.     At  Belvedere  Country  Club,  northeast  of  Hot  Springs, 
the  Bigfork,   tested  by  Albin   (1965)  was   found  to  have  a  transmissi vi ty 
of  2,670  ft/day   (9.42X10'^  m/s).     Many  of  the  springs  in  the  area  issue 
from  the  Bigfork  and  many  of  the  cold-spring  emergences  are  controlled 
by  contact  of  the  Bigfork  with  adjacent,   less  permeable  formations. 
This   association  of  cold  springs  with  the  Bigfork  was  noted  by  Purdue 
and  Miser  (1923). 

The  Arkansas  Novaculite  is   composed  of  three  members  — an  upper  and 
a  lower  member  of  novaculite,   and  a  middle  member  of  chert.     Locally, 
the  upper  part  of  the  formation   is   composed  of  silt-sized  siliceous 
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particles   and  possesses   intergranular  permeability.     The  lower  member 
is   generally  massive  and  dense,   but  is   \/ery  closely  fractured.     The 
middle  member  is   composed  mostly  of  black  shale  and  thin  chert  beds. 
The  Arkansas  Novaculite  is  not  as  permeable  as   the  Bigfork  Chert,  but 
is   locally   intensely  jointed.     Some  cold  springs   issue  from,  and  many 
water  wells   tap,   the  Arkansas  Novaculite. 

The  Hot  Springs  Sandstone  is  a  massive,  quartzitic  sandstone. 
Fairly  large  joints   and  fractures,   as   in  the  novaculite,  create  some 
highly  permeable  conditions,  such  as   at  the  hot  springs. 

The  hot  springs  emerge  from  the  Hot  Springs  Sandstone  near  the  axis, 
on  the  northwest  limb  of  a  southwestward  plunging  anticline.     The 
springs  emerge  between   the  traces   of  two  thrust  faults,  parallel    to  the 
axis  of  the  anticline   (fig.   3).     The   locations  of  the  hot  springs,  shown 
in  figure  2,  generally  lie  along  several   northeast  trending  lines. 
According  to  Bryan  (1924),   these   lineaments  were  inferred  to  be  the 
traces  of  fissures  by  R.    R.   Stevens,  who  first  noted  their  alinement  in 
1890.     Jointing  is   common  in   the  few  exposures  of  the  Hot  Springs  Sand- 
stone in   the  hot-springs-discharge  area.     Thus,   the  hot  springs   are 
associated  with  thrust  faults,  and  normal    faults  and  joints  on  the 
plunging  crestline  of  the  anticline.     Upward  movement  of  the  hot  waters 
from  depth  is  probably  along  the  permeable  fault  zones.     These  fault 
zones  probably  carry  water  to  near  the  surface,  where  the  water  follows 
permeable  joints   to  the  spring  outlets. 

Geologic  sections   in  figure  3  show  the  geologic  structure  in  the 
vicinity  of  the  hot  springs. 
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CHARACTER  OF  THE   SPRING  AND  WELL  WATERS   IN   THE  HOT-SPRINGS  AREA 

Physical   Quality 
Flow  of  the  Hot  Springs 

Flow  of  the  hot  springs  has   been  measured  at  infrequent  intervals 
by  several    investigators.     The   first  to  attempt  measurements  was 
William  Dunbar,   of  the  Dunbar  and  Hunter  expedition,   in   1803,  who 
measured  only   the   largest  springs.     In   1860,  Glasgow  determined  the  hot- 
springs   flow  to  be  450,480  gpd  (1.98  X  10"^  mVs).     Probably  the  first 
accurate  measurement  of  the  total   sprina  flnw  was  made  by  W.  H.  Weed, 
in   1902.     Weed,   after  measuring  or  estii.iating  the   flow  of  each  spring, 
found  the   total    flow  to  be  850,000  gpd  (3.73  X  10"^  mVs).     The  flow  of 
all   but  a  few  of  the  hot  springs   is  piped  into  a  central-collecting 
reservoir.     Hamilton,   in   1931,  calculated  the  maximum  rate  of  filling 
of  the  central-collecting  reservoir  to  be  960,000  gpd  (4.21   X10"2  mVs). 
Park  Superintendent  Libbey,   in  1945,   recorded  that  the  central-collecting 
reservoir  filled  with   15   feet  (4.57  m)  of  water  in   7  hours,  which 
represents   an  average   flow  of  800,000  gpd  (3.51   X   10"^  m^/s). 

The   volume  of  the   reservoir  per  unit  change   in  depth  is   known  and 
provides   a  means   for  periodic  calculation  of  spring  flow.     One  large 
spring  emerges   at  the  bottom  of  the   reservoir.     The  flow  of  this  spring 
decreases   as   the  stage  in   the   reservoir  rises.     Also,  above  a  stage  of 
16  feet,  overflow  occurs.     Flow  calculations   are  thus  made  during  the 
filling  cycle  below  a  stage  of  16  feet  (4.88  m).      Further,   to  avoid 
variable  effect  of  stage  in   the  reservoir  on  the  flow  of  springs  in  the 
reservoir,   the  flow  calculations   are  adjusted  to  a  stage  of  15  feet 
(4.57  m)   by  means  of  a  rating  curve. 
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A  hydrograph  of  the  average  of  calculations  of  monthly  flow  since 
September  1970  is  shown   in  figure  4.     Spring  discharge  has   ranged  from 
750,000  gpd   (3.29  X   10-2  ^3/^)    ^q  q^q^^^  950,000  gpd   (4.17  X   lO'^  mVs). 
Fluctuations  of  spring  discharge   follow  a  seasonal   cycle  during  the  year; 
discharge  is  high  in  the  winter  and  spring  and  is   low  in  the  summer  and 
fall.     One  would  expect  the  flow  of  an  individual   spring  to  fluctuate 
seasonally,   as   does   the  total   spring  flow.     The  flow  of  the  lower  display 
spring   (number  32)   is   shown  by  hydrograph  in   figure  4.  ^ 

In  addition   to  seasonal    changes  in  flow,   individual   springs  show 
long-term  changes   in   flow.     Some  springs  have  ceased  to  flow  and  new 
springs  have  emerged  during  the  past  170  years.     These  long-term  changes 
in  spring  flow  have  not  been  systematically  documented;  many  changes  have 
doubtlessly  passed,   unrecorded.     Measurements   of  Weed  (1902),  Hamilton 
(1932),  and  the  present  study  afford  samples   for  comparisons  of  long- 
term  variations. 

Many  changes  in   individual   spring  flow  occurred  between  1901   and 
1931    because  of  excavation  and  construction  at  spring  outlets,  opening 
new  springs,   and  drilling  the   Fordyce  well    (fig.   2).     Hamilton  noted 
that  the  flow  of  springs   in  group  1,  which   includes   those  at  higher 
elevations,  declined  in  flow  from  168,000  gpd  (7.37  X  10"^  mVs),   in 
1901,    to  124,000  gpd   (5.44  X   lO'^  mVs),   in   1931    (table  2). 

Another  group  of  springs,  at  lower  elevations   along  the  base  of 
a  tufa  cliff,  showed  practically  no  change  in  flow  from  1901    to  1931. 
These  springs   (numbers   2,   4,   6,  9,   11,   13,   15,   17,   19,   and  two  unnumbered 
springs)    discharged  313,500  gpd  (1.376  X  10"^  mVs)   in   1901   and  315,000 
gpd  (1.382  X   lO'^   mVs)    in   1931. 
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Figure  4. —Temperatures  and  flow  of  the  hot  springs. 
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Table  2. — Flow  of  the  hot  springs  in  group   1  in  1901  and  1921 

lAfter  Hamilton,  1932] 


Flow,  in  gallons  per  day 


Spring  number 


1901 


1931 


1 

3,  5,  and  8 

7 r 

10 

22 

23  and  24 

26 -- 

27,  28,  and  29 

47  and  48 

49 

Maurice  Spring 

New 

Total 

■"■  Not  measured  individually. 


28,800 

9,600 

39,218 

21,800 

18,516 

1,760 

18,514 

14,400 

1,723 

2,460 

10,800 

5,000 

25,847 

10,950 

24,418 

{') 

Nonexistent 

13,500 

Nonexistent 

(M 

Nonexistent 

{') 

Nonexistent 

2,400 

167,836 

124,000 
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Thirteen  springs   that  were  measured  by  Weed  in   1901  were  not 
included  in  the   1931    collection  system  described  by  Hamilton   (1932). 
Two  of  these  springs  are  on  the  Arkansas  Rehabilitation  Center,  one  of 
which,  number  39,   is   used  for  supplying  hot  water  to  the  Rehabilitation 
Center.     Four  springs  measured  in   1901  were  dry  in  1931.     Locations  of 
six  springs  were  unknown  in   1931.     Presumably,   these  six  springs  were 
dry   or  had   insignificant  flow   rates   in   1931. 

Hamilton   lists   five  springs    in  his   1931   collection  system  that  were 
nonexistent  in   1901.     Two  of  these   (47  and  48)   have  declined  in  flow 
since   1931   and  one   (49)    is  now  one  of  the  larger  springs   in  the  system. 

The  collection  system  has  not  changed  since  1931,   and  no  springs 
have  been  excavated  nor  have   hot-water  v;g11s   been   drilled,   resulting  in 
a  relatively  stable  period  of  spring  locations.     However,   changes   in  flow 
rates  since   1931    have  been  noted. 

Flows  of  several   springs  were  measured  or  estimated  in  1972  while 
the  springs  were  uncovered  for  sampling  in  January   1972.     These  measure- 
ments showed  a  general   decline  of  flow   in   the  springs   at  higher  elevations. 
Spring  number  1   had  a  very  small    or  no  flow— a  decline  from  28,800  gpd 
(1.26X10-^  mVs),   in   1901,   to  9,600  gpd   (4.21X10"'^  mVs),   in   1931.     Spring 
numbers   47  and  48  declined  from  13,500  gpd   (5.91X10"^  mVs),   in   1931,   to 
8,600  gpd  (3.77X10"^  mVs),   in   1972.     However,   spring  number  7  declined 
from  18,516  gpd  (8.12X10"^  mVs),   in    1901,   to  1,760  gpd  (7.72X10"^  mVs), 
in   1931,  but  rose   to  2,800  gpd   (1.27X10"'^  mVs),   in   1972.     Spring  number 
49,  which  was  nonexistent  in   1901,   had  a  flow  of  58,000  gpd  (2.54X10'3  mVs), 
in  1972. 
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Temperatures   of  the  Hot-Springs  Water 

Temperature  measurements  were  the  first  scientific  data  collected 
on  the  hot  springs.     Dunbar  and  Hunter,   in   1804,   recorded  67.8°  C   (Celsius) 
(154. 0°F)   for  the  hottest  spring   (Weed,   1902).     In   1860,   the  highest 
temperature  measured  by  Owen  was  64.4°C   (148.0°F).     Glasgow  also  measured 
temperatures   in   1860  and  recorded  a  maximum  of  65.6°C  (150.0°F).     The 
maximum  temperature  measured  by  the  Geological   Survey,   in  1972,  was 
61.8°C   (143. 2°F). 

Measurements  of  temperatures  of  individual   springs  by  several   investi- 
gators  from  1890  to  1953  (table  3)   show  maximum  temperatures   of  63.9°C 
(147. 0°F),   in  1901,  64.4°C   (148. 0°F),   in  1931,  and  63.3°C   (146. 0°F),   in 
1952.     Thirteen  of  the  same  hot  springs  were  measured  by  Haywood  in  1901; 
Hamilton,   in  1931;  and  Kuroda,   in  1952.     The  average  temperatures  of  these 
hot  springs  in  1901,   1931,  and  1952  were  58.2°C   (136. 8°F),   57.3°C  (135. 2°F), 
and  58.9°C   (137. 8°F),   respectively. 

These  particular  data  indicate  that  there  has  been  a  slight  decline 
in  maximum  water  temperature  with  time   (0.6°F  in  10  years).     They  are  not 
conclusive  evidence  of  a  general   decline  in  temperature  because  of  dif- 
ferences  in  sampling  points,  variations  in  temperature  with  flow  rates 
of  individual   springs,  and  temperature  fluctuations  in  spring  flow  due 
to  air  temperature. 

Temperature  fluctuations   of  individual   springs  exhibit  an  annual 
cycle  that  seems   to  be  in   response  to  the  annual   cycle  in  air  temperature. 
Graphs   of  the   temperatures   of  four  springs  are  shown  in  figure  4.      In 
some  springs  short-term  fluctuations  superposed  on  the  annual   cycle  are 
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Table  3. — Temperatures ^    in  degrees  Celsius^   of  the  hot  springs 


Spring 


Date  of  measurement 


1890 
(from 
Branner, 
1892) 


1900 
( from 
Haywood, 
1902) 


1901 
(from 
Haywood, 
1902) 


1931 
(from 
Hami  Iton , 
1932) 


1952 
(from 
Kuroda, 

1953) 


1972 

(present 

study) 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


62.6 


59.6 


59.6 
63.6 


46.6 


61.6 


61.9 
51.9 
61.7 
55.9 
61.4 
57.5 
60.1 
35.2 
61.1 
57.4 

36.2 

60.9 
63.9 
60.8 
55.4 
57.3 
56.4 
46.3 
43.3 
57.1 
62.0 
62.3 
62.7 
63.4 


57.1 
53.9 
51.4 
46.0 
48.3 
47.9 
39.0 
48.9 
52.9 


61.7 

53.9 

61.3 

52.4 

61.9 

58.3 

60.8 

36.2 

62.4 

57.2 

56.8 

36.2 

56.3 

62.8 

63.9 

60.9 

56.4 

57.3 

56.1 

44.5 

46.0 

56.5 

62.4 

60.3 

62.9 

61.4 

51.9 

59.8 

57.8 

51.9 

51.4 

46.5 

49.2 

47.3 

43.0 

48.8 

52.6 


54.5 
60.0 
57.2 
61.1 
57.2 

61.1 
62.2 
60.0 
64.4 


64.4 

55.6 
52.8 

44.7 

56.1 
50.0 
57.8 

63.3 


54.4 
45.0 


61.2 

61.2 

61.2 

56.9 
59.6 


59.6 
54.3 
63.3 

59.2 
60.0 
61.1 

54.0 
52.5 

57.5 


62.0 

54.6 

62.2 

52.2 

61.7 
58.6 
62.2 

59.3 

61.2 

56.0 


56.2 


52.5 
57.6 
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Table 

3.—Tempera-i 

tuvesjin  degrees   Celsius ^of  the   hot  springs-' 

-Continued 

Date  of  measurement 

Spring 

1890 
(from 
Branner, 
1892) 

1900 
(from 
Haywood, 
1902) 

1901 
(from 
Haywood, 
1902) 

1931 
(from 
Kami Iton, 
1932) 

1952 
(from 
Kuroda, 

1953) 

1972 

(present 

study) 

38 

58.8 
61.4 
48.9 
46.8 

59.8 
61.4 
48.9 
48.3 
58.3 

60.3 

39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
Maurice 

51.6 

60.6 
50.0 

60.8 

61.3 

46.1 

8.0 

13.0 

51.5 

57.2 
58.6 

58.3 

61.7 
62.2 
61.0 

60.0 
61.8 

58.9 
62.8 

53.3 
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evident  from  detailed  temperature  records.     The  short-time  fluctuations 
are  due  to  external   environmental    factors.     The  Fordyce  Spring  (Spring 
No.   46,   fig.   2)   temperature,   for  example,   responds    to  such  factors  as 
heat- load  changes,   caused  by  opening  the  air-tight  seals  on   the  chamber 
enclosing  the  spring,   and  effects   of  mixing  of  the  hot-springs  water 
with  seepage   from  nearby   rainfall. 

It  v/ould  be  expected  that  spring  temperature  increases  with  increase 
in  spring  flow.     This    relation  has  been  examined  by  scatter  diagrams 
drawn  between   flow  and  temperature   for  the  display  springs.     A  direct 
relationship  is  obscured  because   the  effect  of  ambient  air  temperature 
varies  seasonally  and  the  high  seasonal   air  temperature  occurs  during 
the  period  of  low  spring   flow. 

Silica  Concentration  as  a  Temperature  Indicator 

The  solubilities  of  si  lea  minerals   increase  with  increase  in  tempera- 
ture.    Silica  goes   into  solution   to  the  solution  saturation  points 
hundreds  of  times   as   fast  as   it  precipitates   from  oversaturated  solutions. 
Thus,  a  thermal  water  will    dissolve  silica  minerals   as   its   temperature 
rises,   but  as   the  water  cools   it  will   not  rapidly  lose  silica.     The 
silica  concentration  of  water  can  be  used  as   a  measure  of  the  maximum 
temperature   reached  by   the  water  (Fournier  and  Rowe,    1966). 

The    actual    silica  concentration  of  water  depends  on  the  particular  silica 
mineral    to  which  the  water  is  exposed,   as  well   as   to  the  temperature. 
In   the  hot-bprings  system,   the  main  sources   of  silica  are  the  Bigfork 
Chert  and  the  Arkansas  Novaculite.     Chert  and  novaculite  are  both  composed 
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of  chalcedony,  a  microcrystalline  variety  of  quartz  that  has  its  own 
regular  chemical  and  physical  properties.  The  solubilities  of  chalcedony 
and  quartz  (R.  0.  Fournier,  oral  commun.,  1972)  are  shown  in  figure  5. 

The  analyzed  silica  concentrations  of  the  samples  collected  for  this 
study  are  also  shown  in  figure  5.  The  cold  springs  and  wells  are 
oversaturated  with  respect  to  quartz,  but  are  undersaturated  with 
respect  to  chalcedony.  Though  it  is  likely  that  complex  silica-bearing 
clays  or  other  minerals  influence  the  silica  concentration  of  these  waters, 
it  is  also  possible. that, because  of  their  low  temperatures,  these  waters 
have  not  had  sufficient  time  to  reach  saturation  with  the  chalcedony. 

The  hot-springs  waters  have  a  uniform  silica  concentration  (41-42 
mg/1).  They  are  oversaturated  with  respect  to  chalcedony  at  their  tempera- 
ture when  collected,  but  the  higher  temperature  springs  approach  saturation. 
The  constant  silica  concentration  of  the  springs,  together  with  the  similar- 
ities in  their  other  chemical  properties,  strongly  suggests  that  all  the 
hot  springs  emanate  from  a  single  source  having  a  uniform  geothermal 
environment.  From  the  silica  concentration  (41.5^0.5  mg/1),  the  tempera- 
ture of  the  source  can  be  calculated  as  63.2+0.5°C  (145.8±0.9°F).  The 
different  temperatures  of  the  individual  springs  are  presumably  a  result 
of  differing  flow  paths  and  rates  of  approach  to  the  surface  from  their 
common  source. 

Silica  analyses  were  made  on  groups  of  samples  collected  in  1901, 
1952,  and  for  the  present  study  in  1972.  These  analyses  provide  a  firmer 
base  from  which  to  judge  temperature  trends  with  time.  The  silica  analyses 
of  springs  in  a  given  group  show  less  variation  than  overall  spring 
temperatures.  As  an  indicator  of  maximum  temperature,  the  silica-computed 
temperature  is  not  affected  by  external  factors  such  as  air  temperature 
and  flow  of  the  springs, 

26 


q; 


ex: 

UJ 

tn 

<c 


70 

0 

1 

50 

- 

40 

- 

30 

- 

20 

- 

TEMPERATURE,  IN  DEGREES  CELSIUS 

15  30  45  60 


15 


«3: 
on 


o 


o 
o 


CM 

o 


oo 


10    - 


5    - 


3.7 


3.5  3.3  3.1  2.9 

l.OOO/TEMPERATURE,  IN  DEGREES  KELVIN 


2.7 


Figure  5.— Variation  of  dissolved-silica  concentration  with  temperature, 


Silica  concentrations  and  maximum  temperatures  calculated  from  the  silica 
concentrations   (Fournier  and  Rowe,   1966),  as  well   as  maximum  temperatures 
recorded  at  various   times,  are  shown  in   the  following  table. 


Year 

Maximum 
temperature 
measured 
(°C) 

Average 

Si02 
concentration 

(mg/1) 

Number 

of 
springs 
sampled 

Temperature 
calculated 
from  Si02 
concentration 

1804 

67.8 



— 



1860 

64.4 



— 



1860 

65.6 



— 



1901 

63.9 

46.6 

40 

68.5 

1931 

64.4 



-- 



1952 

63.3 

42.7 

11 

64.3 

1972 

61.8 

41.7 

9 

63.2 

These   results   are  also  plotted  in   figure  6.     According  to  the 
silica  data,   the  spring-source  temperature  has   been  decreasing  at  about 
0.077°C   (0.14°F)   per  year  since   1901.      Further,  measured-maximum 
temperatures   since   1931  — a  period  in  which   the  springs  were  undisturbed- 
decreased  at  about  the  same   rate   (0.073°C,   or  0.13"F  per  year).     This 
coincidence  may  be   a   fortuitous   circumstance  of  sample  timing  and  dis- 
tribution.     In  addition  to  reflecting  the  change   in  the  source  tempera- 
ture,  variations    in  measured  temperatures   reflect  the  sampling  distribution, 
seasonal   changes   in  temperature,   and  changes   in  temperature  due  to  changes 
in   flow  distribution— both  natural   changes   in   flow  distribution   and  those 
due   to  man's  efforts   to  enclose  and  control    the  springs.     Thus,   although 
the  temperature  observations  may  be  subject  to  sampling  bias,   the  silica- 
concentration  trend  of  the  hot  springs   indicates  a  decline  in  temperature. 
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Temperatures  of  the  Cold  Springs  and  Wells 

The   temperatures  of  ground  water  in  wells  in   the  vicinity  of  Hot 
Springs,   other  than  the  hot  springs   themselves,   range  from  12.0°C  to 
52.8°C   (tables  4,5;  fig.    7).     The  warmer  of  these  temperatures  are  the 
temperatures   of  ground  waters   in   the  immediate  vicinity  of  the  hot 
springs.     The  well   on  the  grounds  of  the  Arkansas  Rehabilitation  Center 
has   a   temperature  of  52.8°C   (127. 0°F);   the  water  from  a  well   on   the 
Arlington  Hotel   lawn,  used  to  supply  water  to  the  cooling  plant,  has   a 
temperature  of  32.2°C  (90.0°F);  the  temperature  of  the  well    at  the 
Arlington  Hotel    has   a   reported  temperature  of  22.8°C   (73.0°F).     The  high 
temperatures  are  associated  with  the  abnormally  high  geothermal   gradient 
caused  by  heat  conveyed  to  the  surface  by  the  waters  of  the  hot  springs. 

The  temperatures  of  ground  waters  in  and  near  areas  considered  as 
potential   for  recharge  to  the  hot-springs  system  (outcrop  area  of  Big- 
fork  Chert  and  Arkansas  Novaculite)    range   from  15.0°C   (59.0°F)    to 
26.8°C   (80.2°F).     These  temperatures,   considered  in  relation  to  the 
average  annual   air  temperature  of  17.7°C  (63.9°F),  indicate  that  some 
of  the  waters   in   and  near  the  recharge  area  are  heated  by  geothermal 
processes.      Ground-water  temperatures,   in   the   range  of  15.0°C   (59.0°F) 
to  16.7°C   (62.1°F),   indicate  a  minimum  of  geothermal   heating  and  reflect 
the   fact  that  most  recharge  occurs  during  the  cooler  seasons  of  the 
year.     These   temperatures   indicate  a  rapid  rate  of  movement  of  water 
downward  from  the  surface  and  proximity  to  the  recharge-source  area. 
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Table  5. — Temperatures   of  water  in  the  welts 
[Thermal  gradient  =  (temp.  (°C)-15)/depth  of  well  (ft)] 


Well    no. 
(fig.    7) 


Temper- 
ature 
(°C) 


Date 
of 
measure- 
ment 


Depth 

of 

well 

(feet) 


Thermal 

gradient 

(°C/foot) 


Geologic  formation 


Wl- 
W2- 

W3- 
W4- 
W5- 
W6-- 
W7- 
W8- 
W9- 
WIO- 
Wll. 
W12- 
W13- 
W14- 
W15- 
W16- 
W17- 
W18- 
W19- 
W20- 
W21- 
W22- 
W23- 
W24- 
W25. 
W26- 
W27. 
W28. 
W29' 
W30' 
W31 


15.6 

17.8 

16.4 

16.1 

16.9 

16.7 

15.6 

16.1 

15.3 

16.7 

16.7 

17.4 

15.3 

15.6 

17.6 

15.6 

15.0 

15.8 

16.1 

16.1 

20.4 

16.7 

16.1 

18.8 

18.6 

18.0 

22.8 

52.8 

32.2 

18.0 

16.7 


6-29-72 
9-13-72 
6-  6-72 
7-14-72 
7-14-72 
7-14-72 
6-  6-72 
6-  6-72 
6-  7-72 
6-  7-72 
6-  6-72 
9-25-72 
9-12-72 
6-  7-72 
9-25-72 
6-  2-72 
6-  2-72 
6-  6-72 
6-  2-72 
6-  2-72 
9-27-72 
11-  2-71 
11-  2-71 
1-22-72 
1-28-72 
1-24-72 
3-  8-70 
8-  8-72 
3-  4-70 
9-22-72 
9-22-72 


26 
30 
61 

170 
92 
30 
89 
93 
42 
20 

101 
44 
67 

126 

140 


38 

78 

235 

263 

120 

90 

46 

300 

89 

28 

202 

336 

200 

300 

140 


0.023 
.093 
.023 
.065 
.021 
.057 
.0067 
.012 
.0071 
.085 
.017 
.054 
.0045 
.0048 
.019 


0.0 
.010 
.0047 
.0042 
.045 
.019 
.024 
.063 
.040 
.11 
.039 
.11 
.086 
.010 
.012 


Bigfork  Chert 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 

Do. 
Arkansas  Novaculite 

Do. 
Bigfork  Chert 

Do. 

Do. 

Do. 

Do. 
Hot  Springs  Sandstone 

Do. 

Do. 
Arkansas  Novaculite 
Bigfork  Chert 
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Low  ground-water  temperatures,   below  the  average  air  temperature, 
occur  in   the  area  northeast  of  the  springs   and  in  the  outcrop  of  the 
Bigfork  Chert  and  the  Arkansas  Novaculite.     Water  temperatures   from 
several  wells  in   the  area  are  higher  than   the  average  air  temperature. 
These  temperatures  indicate  qeothermal   heating  due  to  a   long  residence 
time  of  the  water  in   the  aquifer.     Of  the   five  wells  sampled  having 
temperatures  higher  than  the  average  air  temperature,  four  are  flowing 
artesian  wells. 

The  cold-water  springs   in   the   ared   are  generally  wanner  than  the 
well   waters.     The  spring   temperatures   range   from   IS.G'-'C   (60.1°F)    to 
26.8°C   (80.2°F).     The  designation  of  "cold-wa-f-pr"   to  these  springs  is 
for  the  purpose  of  distinguishing  them  from  the  hot  springs  of  Hot 
Springs  National   Park.     The  warmer  springs   in  this  cold-water  group 
could  be   classified  correctly  as   thermal    springs. 

Chemical   Quality 

The     hydrology  of  the  hot  springs   is  controlled  by  complex  stratig- 
raphy and  structure.     Ground  waters  contain  a  variety  of  chemical 
species   in   solution.     Knowledge  of  their  chemical    types   and  concentrations 
is   useful   in  interpreting  patterns  of  ground-water  movement,  source  of 
the  water,   and  in  determining  whether  there  has  been  a  change  in  the 
chemical   concentration  of  the  water  through  the  years.     Chemical   analyses 
may  also  hold  clues   to  the  future  life  of  the  hot  springs. 
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There   has  been   little  change  in   the  chemical   quality  of  the  hot- 
springs  water  during  the  period  of  record  from  1890  to  1970.     A  small 
change  in  the  silica  concentration  is  indicated  and  its  significance  has 
been  discussed.     Minor  variations  in  other  constituents  are  consistent 
with  normal    variations  in  most  ground  water. 

As  part  of  this  study,  samples  for  chemical   and  isotopic  analyses 
were   collected  in   1972  from  15  cold  wells   and  springs  in  the  Hot  Springs 
region  and  from  9  of  the  hot  springs  themselves.     Table  6  gives   the 
results  of  the  chemical   analyses  made  for  this  study.     The  concentration 
of  certain  dissolved  constituents  may  change  between  the  time  of  field 
collection  and  the  time  of  laboratory  analysis  owing  to  outgassing 
temperature  changes     and  precipitation  of  solids.     To  insure  that  the 
analyses  would  represent  the  chemistry  of  the  water  in  the  field  as 
closely  as  possible,  certain  analyses  were  made  in   the  field  and  parts 
of  the  samples  were  treated  before  being  sent  to  the   laboratory,   to 
prevent  changes   before  analysis. 

Temperature  was  measured  at  the   collection  site  using  thermometers 
readable  to  0.1°C.       pH,   alkalinity,  and  dissolved  oxygen   (DO)  were 
also  measured  as  part  of  the  sample  collection  procedure,  because  they 
are  liable  to  change  by  exchange  with  atmospheric  CO2  or  O2  gas  when  a 
water  sample  is  shipped  or  stored.     The  precautions  mentioned  by  Barnes 
(1964)  were  observed  in   the  pH  and  alkalinity  measurements.     The  pH 
values  are  considered  accurate  to  10.02  units.     The  dissolved  oxygen 
analytical   procedure  was   that  described  by  Brown,  Skougstad,  and  Fishman 
(1970,  p,   126).     A  separate  bottle  of  each  sample  was  acidified  at  the 
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time  of  collection   and   taken   to  the   laboratory  for  analyses   for 
aluminum  (Al),   iron    (Fe),  manganese   (Mn),   zinc   (Zn),   calcium  (Ca), 
magnesium  (Mg),   and  strontium  (Sr).     Analyses   for  these  seven  ions, 
as  well   as   for  sodium  (Na)   and  potassium  (K)  were  made  by  atomic  absorp- 
tion spectroscopy.      All    determinations  were  made   by  using  methods 
standard  in  U.S.   Goelogical   Survey  laboratories   (Brown,  Skougstad, 
and  Fishman,   1970).      The  aluminum  analyses  of  samples  collected  in 
September  1972  v;ere  made  using  a   larger  aliquot  than  was   used  for  the 
analyses  of  the  samples   collected  earlier.     Thus   the  precision  reported 
in  table  6  is   greater  for  the  samples   collected  in  September  1972. 
The  precision  of  all    analyses   is  estimated   to  be  plus   or  minus  one-half 
of  the  last  reported  digit.     That  is,   a  reported  2.1  mg/1   implies   a 
concentration  of  between  2.05   and  2.15  mg/1.      A  reported  0.00  mg/1   means 
that  the  constituent  was  analyzed  for     but  was  not  found  present  above 
0.005  mg/1. 

The  chemical   analyses  given    in  table  6  are  arranged  approximately 
in  order  of  increasing   total    dissol ved-soli ds  concentration.     U.S.  Geological 
Survey  station  numbers  of  the  sampling  points   are  given  in  table  7. 
Two,  possibly  3,  of  the  analyses  given  in   table  6  are  not  of  particular 
use  in  describing   the  chemical    history  of  the  hot-springs-area  ground 
waters.     Music  Mountain  Spring  issues  from  Stanley  Shale,  and  its 
chemical   characteristics,  particularly  its  silica  and  chloride  concentrations, 
are  unlike   those  of  waters  in   the  hot-springs-flow  system.     The  R.   B. 
Yates' well,   although  in   the  geographic  and  geologic  area  of  interest, 
is  very  shallow  and  the  water  has    •  high  nitrate  (NO3)  concentration  (8.3  mg/1) 
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Table  7. — Station  numbers  of  Dells  and  springs  given  in   table  6 
Station  number  Name  of  v/ell  or  spring 

343110093025301 Happy  Hollow  Spring 

343252091002301 Belvedere  Country  Club  well 

343204093005501 Frank  Thompson's  house  well 

343211093011501 Sleepy  Valley  Spring 

343227093002201 Bill  Sargo's  well 

343327093000401 R.  B.  Yates'  well 

342954093070101 Music  Mountain  Spring 

343033092584901 Cluster  Spring 

343105092572001 McLendon  Mineral  Spring 

343130093020801 Gulpha  Gorge  well 

343347092594201 Elizabeth  Brown  well 

343052093040802 Whittington  Avenue  Spring 

343056093040601 Whittington  Park  well 

343231093012801 Echo  Valley  Spring  (Big  Chalybeate) 

343112C93033601 Diamond  Mineral  Spring  (Lithox) 

343051093031101 Maurice  Hot  Spring 

343057093031301 Hot  Spring  No.  17 

343053093031101 Hot  Spring  No.  23 

343051093030901 Hot  Spring  No.  33  (Upper  Display) 

343049093031101— Hot  Spring  No.  46 

343058093030803 Hot  Spring  No.  48 

343055093031201 80,000-gal  Ion  reservoir 

343047093031001 Mot   'i^.:  N-  4?  (Health  Services) 

343035093031001 Hot  Spring  No.  49 
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and  high   chloride   (Cl=6.4  mg/1).     This   chemistry  and  the   location   of 
the  well   S-uggest  that  the  water  is   influenced  by  man's   activity  and 
does  not  represent  the  natural   conditions  of  interest  in  this  section 
of  the   report.     The  water  from  Bill   Sargo's  well,   too,     has   slightly 
high   chloride  and  nitrate  concentrations   (3.6  and  1.7  mg/1,   respectively) 
and  may  also  not  be  completely  representative  of  natural    conditions. 

The  water  from  the  hot  springs   is   distinctive  chiefly  by  its 
relatively  low  mineral   concentration   (table  6).     The  mineral   concentra- 
tion of  the  water  is   low  probably  because  the  rocks  associated  with  the 
hot  springs  are  made  up  of  only  a  few  minerals,  each  of  which  has  a  low 
solubility.     Most  ground  water  in  Arkansas  contains  from  two  to  three 
times  more  dissolved  minerals   than  the  hot-springs  water. 

Geochemistry  of  Hydrogen  and  Oxygen  Isotopes 

The  elements   hydrogen  and  oxygen  have  several   naturally  occurring 
isotopes.     Hydrogen  has  a  stable  isotope  of  mass  1,  which  is   common 
hydrogen  or  protium  (H  or  ^H),  and  of  mass   2,  which  is   deuterium  (D  or 
^H),       A  radioactive   isotope  of  hydrogen,   tritium  (^H),  also  occurs   in 
the  environment  and  is   discussed  below.     Common  oxygen  has   a  mass  of 
16  (^^^0),  but  oxygen  of  mass   17  (^■^0)   and  mass   18  (^^0)   also  occurs. 

Absolute  isotopic  abundances  or  ratios     are  difficult  to  analyze 
with  precision,  so  it  is   customary  to  measure  and  express  isotopic 
variations   as  deviations  from  an  arbitrary  standard.     These  deviations 

arc  i^\pr(^ssoil    in   deUd  notation  where 

/r  samp  It: 
^       \  R  standard       , 
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Here,  R  is  the  isotopic  ratio  (D/H.  ^^O/^^O)  and  5  is  in  parts  per 
thousand,  or  per  mil  (%.  )  The  delta  notation  is  less  difficult  to 
use  than  it  might  seem. 

The  standard  to  which  natural-v/ater  isotopic  measurements  are 
referred  is  Standard  Mean  Ocean  Water  (SMOW)    (Craig,   1961a).     The  con- 
centrations of  the  various   isotopic  molecular  species  in  this  standard 
are: 

Ho^sq   s  2,000  ppm,  H.,1    0  =  420  ppm,   and  HD^^^O  =   316  ppm. 
H^^'^O  makes   up  the  remainder  and  all    other  species   total    less   than  one 
part  per  million   (ppm).      Knowledge  of  the    -^0  concentration  of  a  sample 
provides  no  more  information   than   knowledge  of  the   ^^0  concentration; 
therefore,   the    ^'0  concentration   is   generally  not  reported.     The  stable 
isotopic  chemistry  of  water  is   thus   expressed  in   terms   of  60  and  6^^0. 

Natural  waters  are  generally  depleted  in  D  and   '^0  relative  to 
SMOW— that  is,   their  :D  and  ,'^^0  values   are  negative.     This  depletion 
occurs  because   the  vapor  pressure  of  water  molecules  containing  the 
heavier  isotopes   is   slightly   less   than  that  of  common  water,  H2^^0. 
During  evaporation  and  condensation  in   the  hydrologic  cycle,  molecules 
containing  heavier  isotopes   are  concentrated  in  the  liquid  phase.     As 
water  evaporates   from  the  ocean,   the  vapor  is   depleted  in  D  and  ^^0 
and  the  amount  of  depletion  becomes  greater  as   the  temperature  of 
evaporation  decreases.      Further  isotopic  fractionation  takes  place  as 
water  is   condensed  and  reevaporated  during  atmospheric  transport,  and 
the  amount  of  fractionation  is   inversely  proportional    to  temperature. 
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The  D  and   ^^0  concentration  of  meteoric  water— that  is,  water  of 
recent  atmospheric  origin  — varies   regularly  over  the  land  surface  of 
the  earth.     A  great  number  of  measurements  show  that  for  meteoric  water 
not  subject  to  much  evaporation,   cD  and  6^^0  are   related  by  the  expression: 

6D=8c,i"0  +  10.        (Craig,   1961b)  (1) 

Further,   the  amount  cf  depletion   increases  with  altitude,   latitu^p,  anH 
distance  from  the  ocean.     The  general   pattern  of  isotopic  distribution 
in  North  America,   at  least,  has  been  napped   (Sheppard,  Nielsen,   and 
Taylor,    1969). 

Figure  8  shows    the   variation   of  >.'D  and  6  ^  ^0  from  several   hydro- 
thermal    localities   and  the  trend   line   for  meteoric  waters,  equation   1. 
The  meteoric  waters   from  the   localities  shown  illustrate  the  tendency 
for  increasing  depletion  from  the  oceanic  isotopic  composition  with 
latitude,   altitude,   and  distance. 

Figure  8  illustrates   that  in  many  hot  ground-water  systems   there 
is  enrichment  in   the   ^^0   concentration  of  the  waters.     An  explanation 
for  this  shift  might  be   that   the  hot  waters   represent  mixing  between 
meteoric  waters  and  waters   emanating  from  hot  or  molten   rock  at  great 
depths  within   the  earth — that  is,  magmatic  water,   the  isotopic  composi- 
tion  of  which   is   not  precisely  known,   but  which  may  well   be   in   the  range 
shown   in   the  figure.     Were  this  explanation  correct,   there  should  be  a 
shift  in   5D,   as  well    as   in  6  ^  ^0  toward  whatever  isotopic  composition 
may  be  typical   of  magmatic  waters.     As  is  particularly  well   illustrated 
by  the  Nevada  and  Yellowstone  Park  series,  such  a  deuterium  shift  does 
not  occur.     Because  a  deuterium  shift  is   absent  and  because  of  the 
existence  of  another  process    to  explain   the   ^^0  shift,  magmatic  waters 
are  probably  not  significant  contributors   to  hydrothermal   systems. 
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The  process   bringing  about  the    ^^0  shift  is  probably  an   isotopic- 
exchange   reaction  between  the  water  and  rock  in  a  system.     Such  a 
reaction     for  a  silica  mineral     can  be  written: 

S-j]GoiJ-.(;,i-Hoi^'0>Siif02+H2^^0. 
Because    ^ '^O   in  most  igneous  minerals   is  enriched   relative  to  SMOW,   the 
effect  of  this   reaction   is   to  drive  the  6^^0  in  exchanging  waters   toward 
positive   values.     Further,  because  of  the  low  hydrogen  concentration  of 
most  of  the  rock,   the  6D  of  waters   is   relatively  unchanged, 

Isotopi c-exchange  reactions,  as  in  all   chemical    reactions,  proceed 
faster  at  high   temperatures   and   their  effects   are  generally  seen   only 
in  systems   in  which   the  water  temperature  approaches  or  is  higher  than 
100°C   (212°F).     Because  of  the  utility  of  stable  isotopes  in  previous 
interpretative  work  in  such  systems,   D  and   ^^0  measurements  were  made 
on  several   samples   from  the  hot-springs  system. 

Isotopic  analytical    results   are  given  in  table  8  and  shown  in 
figure  8.     Deuterium  in   the  cold  waters   has  an  average  6D  =  -28.6%c  ; 
whereas,   the  average  of  the  hot  springs   is   insignificantly  different  at 
-28.3/0..     These  values  agree  with  previous   estimates   (Sheppard,  Nielsen, 
and  Taylor,   1969)    that  deuterium  in  meteoric  water  in   the  central  Arkansas 
region  should  be  within   a  few  per  mil   of  -30.     The   coincidence  of  values 
for  the  hot-   and  cold-water  sampling  points  is  strong  evidence  that  there 
is  no  detectable  source  of  meteoric  water  contributed  to  the  hot  springs 
from  outside  the  region. 

The  average  6^^0  of  the  hot  springs   is  -5.2%*,  and  the  average  of 
the  cold  springs   is  -5.1%^.      The  lack  of  oxygen-isotope  shift  reinforces 
the  deuterium  evidence   that   there  is   only  a  negligible  amount,   if  any, 
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-■*'  of'e''  than  meteoric  water  in   the  hot  springs.     The   lack   of  oxygen 
sni    t  also  snows   that  t'-.ere  is  no  oxygen-isotope  exchange  between  the 
»"0t     ;n'-;g5   and  oxygen-bearing  minerals  in   the  aquifer.     This   lack  of 
oxygen-i sotooe  exchange  could  be   true  only   if   (1)    the  maximum  tempera- 
ture wiinin   the  hot-springs  system  is   not  much  higher  than   the  springs' 
s-rface  zerrteratures ,   and   (or)    (2)    the  flow  rate  of  the  waters   through 
the.  hottest  Dart  of  the  system  is    relatively  high,   so  there  is   little 
■^.ime   for  excnange. 

In  summary,   the  deuterium  and    '^0  data  suggest  that  the  concept 
given  elsewhere  in   this    report  concerning   the  hydrology  of  the  Arkansas 
hot  springs   is   reasonable  and  that  (1)   virtually  all    the  water  in   the 
hot  springs    is    local,  meteoric  water.      Thus   any   theories   implying  long- 
distance  (t^^at  is,  many   tens   or  hundreds   of  miles)   water  movement,   or 
the  presence  of  juvenile  or  magmatic  water  in   the  springs  must  be  dis- 
regarded; and  (2)   the  maximum  temperature  reached  by  most  of  the  hot- 
springs  water  cannot  be  many  tens   of  degrees  higher  than  the  springs' 
surface  temperature,  and  the   resident     time  of  the  hot-springs  waters 
in   the  heated  part  of  the  system  is   a   relatively  short  time,  that  is, 
not  more  than  perhaps   a  few  hundred  years.     The  silica  concentration  of 
the  hot-sDrngs  water,   discussed  in  a  previous  section,  and  the  ^'^C 
".oncentration,  discussed  below,   refine  these  conclusions   further, 

~he  isotope  of  hydrogen  is   tritium  (H^  or  T),  which  has   an  atomic 
.riass  of  3  and  is    radioactive,  with  a  half-life  of  about  12^^  years. 
Tritium  is   formed  continuously  by  cosmic  rays  impinging  on  the  upper 
atmosphere.     This  natural    tritium  is   present  in  precipitation  at  levels 
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o""  ~rom  1    to  10  T  atoms   per  10^'^  H   atoms.      Tritium  is  measured  by 
ara"!    zinc  -' ts   rate  of  radioactive  decay  in  a  water  sample  and  the  results 
are  expressed  as   tritium  units    (TU),  one  of  which  is  equal    to  a  T/H 
ratio  of  "      ".     Thus,  natural    tritium  is  present  in  the  range  of  from 
1    to  10  TU 

Large  quantities  of  tritium  are  produced  by  nuclear  devices,   and 
the  atmospheric  testing  of  such  devices  — particularly  fusion  devices 
(hydrogen  bombs)  — from  the  early  1950's   through   1962  raised  the  level 
of  tritium  in  precipitation   to  many  times  its  natural    level   of  from  1 
to  10  TU.      Peak  tritium  levels  occurred  in   the  spring  of  1963,  when 
precipitation   at  St.   Louis,   for  example,   reached  levels  of  more  than 
2,500  TU.     Since  then,   tritium  levels  have  been  decreasing  at  about  30 
percent  per  year.     During  the  1950' s,   tritium  levels  were  in  the  range 
of  several    hundred  tritium  units,   also  well   above  natural    levels. 

The  tritium  concentrations  of  samples   collected  from  the  hot  springs 
are  given   in   table  8. 

Results  of  radioisotope  measurements  of  tritium  and  also  of  ^'^C, 
that  are  discussed  below,  are  expressed  with  a  statistical-error  term 
corresponding  to  one  standard  deviation   (la).     There  is  a  probability 
of  2    ,n   3  that  the  true  value  of  a  quantity  is  within  the  lo  range.     For 
Diamond  Mineral   Spring   (table  8),  with  tritium  shown  as  l.StO.A  TU, 
there  is   a  67-percent  probability  that  the  true  tritium  concentration  is 
between  1.1   and  1.9  TU,   and  a  95-percenL  (2a)    chance  that,  it  is  between 
0.7  and  2.3  TU.  =* 

Of  the  nine  hot  springs,   four  have  tritium  levels   that  are  so  low 
that  they  nay,  within  a  95-percent-confidence  limit,   contain  no  tritium; 
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that  is,   the  Iz   range  of  their  analyses  includes  zero.     The  rest  of 
"•^e  not  springs   contain   tritium,  but  not  more  than  3  TU.     A  sample 
that  contained  3  TU  could  be  water  that  fell   as   rain  in  the  1940's 
that  had  an  initial    tritium  concentration  of  about  10  TU.     This  explanation 
is   unlikely,   though,   because  evidence  discussed  below  suggests   that 
the  hot-springs  waters  have  ages  of  thousands  of  years,  and  there  is 
another  explanation  for  the  3-TU  level  which  also  accounts   for  the  dif- 
ferent tritium  concentrations  of  the  several   springs. 

It  is  possible  that  some  of  the  hot  springs  include  a  small   pro- 
portion of  water  of  \i^r)i  recent  origin  which  mixes  with  the  main  flow 
of  heated  water  near  the  springs'    outlets.     Such  water  could  contain 
water  from  the  1960's   that  has   tritium  levels  of  several   hundred  or 
more  tritium  units.     An  admixture  of  less  than  1   percent  of  such  water  would 
raise  the  tritium  level    to  that  measured,  and  such  an  amount  would  be  too 
small    to  measurably  affect  any  other  characteristics  of  the  sample. 
Thus,   the  tritium  level    suggests   that  there  is  no  siigiTificant  amount  of 
water  of  less   than  several   decades  old  in  any  of  the  hot  springs. 

Of  the  14  cold  wells  and  springs,  5  have  tritium  levels  higher  than 
25  TU  and  the   rest  have  levels   lower  than  3  TU.     Those  with  little 
tritium,  such  as   the  hot  springs,   contain  no  significant  amount  of 
recently  recharged  water,  whereas   those  with  high  tritium  concentrations 
contain  significant  quantities  of  water  recharged  within   the  past  15 
years.       Their  high   tritium  concentration  reinforces   the  preceding  con- 
clusion  that  the  chemistry  of  Belvedere  County  Club  well   and  Frank 
Thompson's   house  well    is   that  of  recharge  to  the  system,  and  that  Bill 
Sargo's  well   and  R.    B.   Yates'   well    could  be  influenced  by  man's  activity. 
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The  fact  tnat  Happy  Hollow  Spring,  although  chemically  identical  with 
the  high-tritium  Belvedere  County  Club  well   and  Frank  Thompson's  house 
well,  contains  only  1.110.6  TU  is  somewhat  puzzling.     A  possible  explana- 
tion  is   that  the  water  in  this   spring  fell   as    rain  at,  say,   5  TU  in  the 
prebomb  era   (pre-1953)   and  has  no  component  of  post-1953  water.     Sleepy 
Valley  Spring  shows  only  slight  effects  of  mineral-water  reactions, 
and   its   high   tritium  concentration   is    therefore  not  unexpected. 

Carbonate  Geochemistry 

The  most  common  inorganic,   carbon-bearing  chemical   species   are  carbon 
dioxide   (CO   ),   a  gas  which   readily  dissolves   in  water,   bicarbonate   (HCO,"), 
the  predominant  carbon  species   in  most  waters   in  nature,  and  carbonate 
(COr~)>   found  in  such  minerals  as  calcite   (CaCO^),   but  also  present  in  low 
concentrations   in  solution.     The  prevalent  species   in  a  given  solution 
depend  on  the  hydrogen-ion   (H^)    concentration  of  the  solution.     The  pH 
of  a  solution   is   a  measure  of   its   H"*"  concentration   (pH  =  -log  H"*").     A 
solution  having  a  low  pH  has   a   relatively  high  H"^  concentration,  or  is 
acidic;  whereas,  one  having  a  high  pH  concentration   represents   an  alkaline 
solution,  one  having  a  relatively  low  H'''  concentration. 

The  interreactions  among  the   three  carbonate  species   are: 

CO. (gas)   +  HO  ->  H^COj,  (2) 

HXO  .  -V  H"^    f  HCOr,  (3) 

and  HCOr  ->  H"^  +  003'^,  (4) 

Water  in   contact  with  a  gas  phase  containing  CO^  will   dissolve  CO^, 

according  to  reaction  2,   in  an  amount  proportional    to  the  partial   pressure 

of  CO.    in   the  gas    (P^   ).     Some  of  the  H^COg  so  formed  will   dissociate, 
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by  reaction  3,  Z')   bicarbonate  and  hydrogen  ion,  which  will  lower  the 

pH  of  the  solution.   If  an  alkaline  substance  is  present,  the  H"^  so 

formed  will  be  removed  from  solution,  permitting  more  HCOj"  to  form, 

Carbonate  minerals  are  common  alkaline  substances  and  are  dissolved  by 

H  by  the  reaction 

CaCOj,   ,.^>+H'^-Ca"-+HCO  '.  (5) 

^(solid)  3 

If  a  water  originally  charged  with  CO,  is  brought  into  contact  with 
a  carbonate  mineral,  the  reaction 


CO  +H,0+CaCO ->Ca^^+2HC0  "  (6) 


will  occur.  This  is  the  most  important  reaction  controlling  the  chem- 
ical character  of  waters  in  the  hot-springs  area,  and  together  with 
the  solution  of  silica,  discussed  previously,  also  controls  the  chemistry 
of  the  hot  springs  themselves. 

In  discussing  the  carbonate  chemistry  of  these  waters,  the  term 
"total  dissolved  carbonate"  is  often  used.  This  term  refers  to  the 
molar  sur  of  the  dissolved  carbonate  species  H  CO  ,  HCO  -,  and  CO  "2. 
This  is  not  an  analyzed  value,  but  is  calculated  from  the  HCOs"  3nd 
pH  values  given  in  table  2,  and  the  eqilibrium-constant  expressions  corre- 
sponding to  reactions  2,  3,  and  4.  These  calcultations  were  made  using 
the  computer  program  WATEQ  (Truesdell  and  Jones,  1974).  The  total 
dissolved-carbonate  concentrations  so  calculated  for  the  waters  sampled 
are  given  in  table  9. 

Some  of  the  wells  and  springs  were  sampled  for  analyses  of  their 
dissolved-gas  concentration.  Figure  9  shows  the  relation  between  the 
analyzed  partial  pressure  of  CO^, (^CO,,)  with  the  ^CO,,  calculated  from  the 
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values  of  pH  and  HCO3'  given  in  table  6.  Some  of  the  analyzed  values 
are  larger  than  the  calculated  values,  but  most  agree  within  25  percent. 
Although  this  may  appear  a  large  error,  it  introduces  very  little  error 
into  the  much  larger  calculated  total  dissolved-carbonate  concentration, 
which  is  the  parameter  of  interest. 

Water  entering  a  ground-water  system  is  exposed  to  gases  from  the 
earth's  atmosphere  and  from  the  gas  phase  in  the  soil  zone.  The  soil 
air  is  the  last  gas  phase  to  which  the  water  is  exposed  before  it  enters 
the  water-saturated  zone,  so  the  composition  of  the  soil  air  controls 
the  type  and  amount  of  gases  initially  dissolved  in  ground  water.  The 
main  difference  between  the  open  atmosphere  and  soil  air  is  that  the 
soil  air  contains  considerably  more  CO^.  In  the  open  atmosphere,  the 
partial  pressure  of  CO,  (CO  J  is  about  3X10"-  atm  (atmospheres).  Water 
in  contact  with  this '^CO^  at  17^C  contains  about  1.3X10"^  millimoles  per 
liter  H2CO3,  or  about  0.8  mg  H2CO3/I.   In  soils,  though,  much  of  the 
CO2  is  produced  by  plant-root  respiration  and  by  decay  of  plant  debris, 
so  the  CO-  in  soil  air  is  one  or  more  orders  of  magnitude  higher  than 
in  the  open  atmosphere. 

Because  of  differences  in  soil  and  plant  types,  the  soil  air  CO^  may 
differ  from  place  to  place  but  should  be  reasonably  constant  within  a 
given  area.  An  analysis  of  the  composition  of  the  gas  dissolved  in  Happy 
Hollow  Spring,  one  of  the  waters  containing  dissolved  material  probably 
'lei  ivrd  mainlv  n  oir,  almosphcriv,  prcv-ipi  tdtion,  showed  its  CO^  to  be 
3X10"^  atm,  equivalent  to  an  H  CO  ,  concentration  of  about  1 .4  millimoles 
per  liter,  or  8/  mg/ I .  Cds  andly^es  were  not  made  for  the  other  two 
most  dilute  samples  — from  Frank  Thompson's  well  and  the  IJcl vedere  Courilry 
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Club  well— but  analyses  of  their  HCO3"  and  H"*"   (pH)   concentrations  were 
made.     From  these  data  and  from  the  known  equilibrium  constants   for  the 
reactions  between   carbonate  species   (1    and  2  above)   the     CO2   of  the  gas 
in  equilibrium  with  them  can  be  calculated  (Truesdell   and  Jones,   1974), 
For  the  water  from  Frank  Thompson's  well,   the     CO2   so  calculated  is 
about  2X10"2.     The  average   total    dissolved  carbonate   (H^COg+HCO^")   in 
these  three  wells  is  1.15  millimoles  per  liter. 

The  chemistry  of  the  remaining  cold  waters  sampled  differs  from 
that  of  atmospheric  precipitation,  as   represented  by  Happy  Hollow  Spring, 
Frank  Thompson's  well,   and  the  Belvedere  Country  Club  well,   primarily 
in  having  higher  alkaline  earth   (Ca"*"^,  Mg"*"^ ,  and  Sr"*""^)   and  bicarbonate 
(HCO^")   concentrations.     Reactions   representing  the  solution  of  mineral 
carbonate  are  written   in   terms  of  an   idealized  CaC03  mineral.     In  nature, 
CaCOg  may  contain  other  elements  of  the  alkaline-earth  group  such  as 
magnesium  (Mg     )   and  strontium  (Sr"*"   ),   replacing  calcium.     These  sub- 
stitutions have  little  effect  on   the  carbonate  in   the  ground  water  and 
so  in  discussing  the  measured  water  chemistry,   the  total   alkaline-earth 
concentration   (Ca+Mg+Sr)   is  used,  whereas  in  discussing  reactions,  only 
the  idealized-Ca"^^   form  is   referred  to. 

Figure   10   is   a  graph   showing   the   concentrations   of  alkaline  earth 
against  total    dissolved-cabonate  concentrations   (H2C03+HC03~  +  C03"^). 
On   this   graph,  the  point  at  zero-alkaline-earth  concentration  is   at  the 
average   total-carbonate  concentration  of  the  three  most  dilute  samples 
(1.15  millimoles  per  liter).      From  this  point,   lines  are  drawn  showing 
the  trend  of  compositions  which  would  result  from  solution  of  carbonates 
by   reactions   5  and  6.     Three  of  the  samples —Sleepy  Valley  Spring, 
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Cluster  Spring,   and  Whittington  Avenue  Spring — fall   near  the   line  of 
reaction  6.      In  these  samples,   the  dissolved  carbonate  in  excess  of 
that  in   the   recharge  water  is   half  from  solution  of  carbonate  minerals 
and  half  from  solution  of  further  soil-zone-derived  CO2.       The  other 
samples   plot  between   the  two   lines,   and  a  greater  proportion  of  their 
excess   carbonate  is   from  solution  of  carbonate  minerals   than  from 
further  CO,   solution. 

To  adjust  ^'^C   (radiocarbon)  measurements  on  ground  waters  so  that 
they  can  be  used  to  calculate  ground-water  ages,   it  is  necessary  to 
know   the  proportion  of  total    dissolved  carbonate  that  has  come  from  solu- 
tion  of  soil -zone  CO2.     This   factor  can   be   derived  from  the  chemical 
information   in   figure  10,  but  before  describing  the  method,   a  discussion 
of  the  isotopic  chemistry  of  carbon  is   appropriate. 

Isotopic  chemistry  of  carbon 

The  element  carbon   (C)   has   three  naturally  occurring  isotopes  with 
masses  of  12,   13,   and  14.     Carbon  of  mass   12   {^^C)  and  13   {^'^C)   are  stable 
isotopes,  whereas   ^'^C  is  radioactive.      In  any  chemical    reaction  such  as 
the  solution  of  CO^  to  form  HCO3",   or  the  formation  of  carbohydrates  from 
CO2  by  plants  during  photosynthesis,   the  light  isotope  (^^C)   reacts  faster 
than  the  heavier  isotope   (^^C).     The  isotope  ratios  of  reactant  and  product 
species  are  different,   and  isotopic  fractionation  occurs.     Carbon-isotope 
ratios  are  measured  relative  to  an  arbitrary  standard,   a  calcite  of  marine 
origin   known  as   PDB  and  expressed   in   the  6-notation,   as  are  hydrogen  and 
oxygen  isotope  ratios. 
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<5  ^  ^C  values   for  samples    from  a  single   type   of  carbon-bearing 
material    are   relatively  constant,  but  one  material   may  have  quite  dif- 
ferent 6 '  ^C  values  from  another.     Marine   limestones,   for  instance,  have 
6^^C  values   close   to  zero   (because   a  marine   carbonate  is   used  as   the 
standard);  whereas,  atmospheric  CO     and  plants   typical   of  the  hot-springs 
region  are  depleted  in    -  -C  and  have  6  ^  ^C  values  of  about  -7%c   and  -24%^  , 
respectively.     These  large   differences  make  5 ' ^C  measurements  valuable 
in  determining  the  sources  of  dissol ved-carbonate  species    (Pearson  and 
Hanshaw,   1970). 

The  naturally  occurring  radioisotope  of  carbon   (^''C)   has   a  half- 
life  of  about  5,730  years.      Natural    ' ''C   is   formed  by  cosmic-ray   reactions 
in   the  upper  atmosphere  and  is   found  at  virtually  a  constant  concentra- 
tion  in   atmospheric  CO   .    ^''C  concentrations   are  measured  and  reported 
as   percentages   relative  to  a  standard  whose    ''C  concentration  is  that 
of  an  ideal   plant,   grown   in  1950,   in   the  absence  of  any  atmospheric 
CO-,   of  industrial    origin.      ■''C   is   produced  as   a  byproduct  of  nuclear 
explosions,   and  in   the   last  25  years   the    ^''C   level    in   the  atmosphere 
has   risen  considerably  above  its  natural    level    (100-percent  modern). 

Plants  using  atmospheric  CO2   assimilate   ^"C.     When  a  plant  dies, 
the  amount  of   '^""C  it  contains   decreases  because  of  radioactive  decay, 
and  the  amount  remaining,    relative   to   that  of  an  ideal    living  plant, 
measures   the  time  that  has  elapsed  since  its   death.     This   is   the  prin- 
ciple on  which   the  radiocarbon-  (^""C)    dating  method  is  based. 

The  equation   for  radioactive   decay  is 


ATtd '-  '^P'""' 
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where     sin       is   the   ratio  of  sample  to  standard  activity,   X  is   the  decay 

^td 
constant  of  the   isotope   (1/8,266  for  ^^*))   and  t  is   the  age  in  years. 

(The  half-life  of  an   isotope  is   the  value  of  t  when     sm/^  ^  .=0.5,  that 

is   t-i/2  =  ln   2/a=5,730  years   for   ^''C).      For  calculation,   the  decay 

equation   is  often  written   using  natural    logarithms   (In)   as: 

t  =   -8,266   In  '     ^ 


\Astd/ 

As   discussed  in   the   foregoing,   the   C0;_   in   the  soil    air  is   a 
product  of  plant   respiration  and  decay.      Because  of  its  plant  origin, 
soil-air  CO     has    the  same   carbon-isotopic  composition   as  plants;   its 
6  ^  ^C  value  is   about  -24/o.-    and  its    ^""C  concentration  is   close   to  100 
percent   (in   the  prebomb  era).      Soil    CO2   dissolved  by  ground  water  will 
also  inherit  its  plant  isotopic  character,   and  if  soil   CO2  were  the  only 
source  of  carbonate  dissolved  in  a  ground  water,   its    ^'^C  concentration 
would  directly   reflect  the   age  of  the   ground  water— that  is,   the  time 
elapsed  since  the  water  left  the  soil    zone. 

Ground  water  generally  contains  mineral    carbonate,  which  is   commonly 
of  marine  origin   and  has   a  6^^C  value  near  0°L  .      Because   of  its   great 
age   relative   to  the  half-life  of  ^"^C,  mineral   carbonate  has  no  ^'^C. 
Solution  of  ^'^C-free  carbonate   lowers  the  measured- ^'^C  concentration  of 
carbonate  dissolved  in  ground  waters  and  gives   the  water  a  falsely  old 
apparent  age.      However,   if  the   relative  amounts  of  soil-derived  and 
mineral    carbonate   are  known,   this  effect  can  be   corrected. 

For  example,   consider  a  water  of  a   ^''C  concentration  of  lOO-percent 
modern,   that  is,  an  age  of  zero  years,  with  a  total    dissolved  carbonate 
of  4.00  millimoles   per  liter,   of  which   3.00  millimoles   per  liter  is   from 
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soil   CO^.   and   1.00  millimoles   per  liter  from  carbonate  minerals.     The 
measured  -^''C   concentration   of   this   water  will   be: 

((3.0X1.QH1.QX0.0))   =  75-percent  modern, 

and  its  apparent  age  will    be: 

t  =  -8,266   In   (75/100)   =   2,380  years 
instead  of  its   real    age,   zero  years. 

By  knowing  the  proportions  of  the  sources   of  carbonate  to  a  ground 
water,   it  is  possible  to  adjust   the  measured- ^"^C  concentration  to  its 
correct  value  and  so  obtain   a  corrected  age.     The  adjustment  factor,   P, 
is   the  ratio  of  soil-air-derived-   (=plant)   carbonate  concentration   to 


(7) 


total    dissol ved-carbonate   concentrations: 

P  ^plant 

^total 

For  this  example, 

P=MO  =  0.75. 
4.00 

The  adjusted   ^^C  concentration  is 

^'C   (adjusted)^'"^   (measured)    , 

P 

and  for  this  example 

I'^C  (adjusted)=  ZLjercein  =   iqo  percent, 

0.75 

corresponding   to  an   age   of  zero,   as   specified. 

Chemical    information  on  ground  waters  can  often  be  used  to  determine 

the  proportion,   P,   of  soil-air-derived  carbonate  to  total   carbonate.     A 

graphical    representation  of  changes   in  water  chemistry,   such  as   that  given   in 

in  figure  10,  can  be  useful.   Consider  Whittington  Avenue  Spring.     The  total 
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carbonate  concentration  of  this  sample  (3.85  millimoles  per  liter)  is 

the  sum  of: 


Carbonate  dissolved  by  atmospheric 
precipitation  in   the  recharge  area, 
here  taken  as   the  average  of  the 
three   very   dilute  samples,   Happy  HoV 
low  Spring,   Frank  Thompson's  well, 
and  the  Belvedere  Country  Club  well 
shown   in   figure   10: 


1.15  millimoles   per  liter. 


Additional    carbonate   from  the 
reactions : 


CaCOa+H'^^Ca^^+HCO 


(5) 


and 


CO„+H.O+CaCO,->Ca     +2HC0 

2^3  3 


(6)     2.70  millimoles   per  liter. 

Note  that  in  reaction  5  all  the  bicarbonate  is  from  the  solution  of 
mineral  carbonate,  and  that  an  equivalent  amount  of  alkaline  earth  (Ca) 
also  appears  in  solution.  In  reaction  6,  half  of  the  bicarbonate  is  of 
mineral  and  half  of  gaseous  (  =  soil  air=plant)  origin.  In  reaction  6, 
twice  as  much  bicarbonate  as  alkaline  earth  is  produced,  but  the  amount 
of  mineral  carbonate  added  still  equals  the  amount  of  alkaline  earth, 
as   in   reaction  5. 

The  Whittington  Avenue  Soring  contains  1.34  millimoles  per  liter 
of  alkaline  earth  in  solution,  and  from  the  precedina  paragraph,  this 
will  be  the  amount  of  mineral  carbonate  dissolved,  as  well.  That  is, 
of  the  2.70  millimoles  per  liter  total  carbonate  present  in  excess  of 
that  in  the  recharge,  1.34  millimoles  per  liter  is  of  mineral  origin 
and  1.36  millimoles  per  liter  (2.70-1.34)   of  plant  origin. 

To  sum  these   considerations   algebraically: 

^total~^recharge"^^minpraT^^gaseous    '  ^   ' 

where  C  is   the  carbonaU:   concentrcition  from  the  subscript  sources. 
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^recharge  ^"^  ^gaseous  ^^^   of  plant  (-soil-air)  origin  and  Cminepai  = 
^alkaline  earth-  Eauation  8  becomes: 

^totar^plant'''^6lkaline  earth, 
or 

Slent'^^tctar^alkaline  earth*  ^^^ 

The   ^^C- adjustment  factor,  equation  7,   Com  ncv  be  solved: 

P=^pl_a_nt  -   ^totc^r^alkaline  earth  =1-  /  ^alkaline  earth  \ 
^total  ^total  \^total    carbonate/ 

For  the  V.'hi  ttington  Avenue  Spring, 

The  stable  carbon-isotopi  c  composition   of  the   total    carbonate  dis- 
solved in  a  ground  water  also   reflects   the  proportions   of  the  source 
carbon.     As   discussed,   the  6^'K  value  for  plants  is  about  -24%.,  whereas 
that  for  marine  nineral    carbonate   is   about  0".'   .     The   5^^C  value  of  a 
sample  will    be   in   proportion   to   the   amounts   and   isotopic  compositions   of 
its   carbonate  sources: 

xi3r  =^i3r  I  i'"lf""''t  I +o?r  /^mineral\ 

^sample"     Slant  c~t]    ^     ^mineral       'VrrT     * 
^         \^total/  \   Ctotal    / 

Because   '  ^ 'C^iincpg]   =  01,  ,   the   last  tern  is   zero  and 

6^30         ,   -6'3c  /SlanA       ^  (11) 

sample  plant     r 

\^total/ 

vvhich,   on   recalling  equation   7  and  substituting  an  average  value 
for  i'"Cpi3j,^,   t^ccomes 

°     ^sample 
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Carbon-isotope  data  and  calculated  total   carbonate  concentrations 
of  samples   in   the  hot  springs   area  are  given   in   table  9.      For  each,   a 
value  for  P,   the  proportion  of  plant  carbonate,  calculated  using 
equation   10,  is  shown.     Also,  6^^C  values,   calculated  from  equation 
12,   are   given   for  comparison  with   the  measured-6  ^  ^C   values. 

Agreement  between  measured- and  calculated-6^ ^C  values  is  excellent, 
suggesting  that  the  geochemical   processes,   described  in  the  preceding 
paragraphs,   adequately  model    the  field  situation.     The  average  value 
of  the   differences   between   the  calculated-   and  measured-6 ^ ^C  values   is 
0.75o   ,   and  the  average  value  of  the  measured-5 ^ ^C  values  themselves 
is    16.71..     The  error  associated  with   the   calculated  values   of  P  can 

be  estimated  as  Q»  '^^  X  100  =  4  percent.     This   is   larger  than  the   ^^C- 

16.7 
analyzed  errors  and  corresponds   to  an  error  band  of  i500  years   for  a 

sample  of  4,000  years  and  of  ±1,000  years   for  a  sample  of  9,000  years. 

The  adjusted  ages  of  the  cold  wells   and  springs  that  were  sampled 
are  consistent  with   inferences  previously  drawn   from  their  chemical 
character  and  their  tritium  concentrations.     The  Belvedere  Country  Club 
well    and  Frank  Thompson's  well,  which   represent  water  recharging  the 
system,   have   ^^C   concentrations   near  100-percent  modern.     The  R.   B. 
Yates'   well,  which  on   chemical    criteria   is  suspected  to  include   recent 
water  showing  man's   influence,  has  both   ^''C  and  tritium  levels   found 
only  in   the  past  20  years. 

The  adjusted  ages  of  the  remaining  cold-water  samples   range  from 
3,700  to  8,900  years  and  increase  in  age  with   increasing  total   carbonate 
concentration. 
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Age  of  the  hot-springs  water 

As   the  analyses   in   table  6  show,   the  hot  springs  have  almost 
identical    chemical   characteristics,   suggestinq  that  they  have  a  common 
origin.     The  only  detectable   differences   among   them  are   in   their 
carbonate  concentrations,  shown  by   the  different  HCO/  and  pH  values 
in  table  5,   and   in   the   total    carbonate  concentrations   calculated  from 
them  and  given   in   table  9. 

Calculations   have  been  made  to  determine  whether  the  waters  of 
the  hot-sprinas   system  are  saturated  vn  th   respect   to   the  mineral    calcite 
(CaCO.)  — that  is,   to  determine  whether  the   calcium  and  carbonate   con- 
centrations of  the  waters   are  such   that  they  will    dissolve  more  calcite 
or  deposit  calcite   fror  solution.     The  mineral    calcite   (CaCO^)   dissolves 
according  to  the   reaction 

CaCO  ,  -Ca"^  ^+C0  "  I 

3 

At  saturation,   the  product  of  the   thermodynamic  activities   of  the 
dissolved  naterials     a^  +2XaQQ  "       equals   a  constant  (KT),   the  value  of 
which  is   a  function   of  temperature.     That  is,   at  saturation 

^Ca"*'2^^C0"'   "  ^^  ' 
It  is  convenient  to  express   the  degree  of  saturation  of  a  water  in  terms 
of  the   logarithm  of  the   ratio  of  the   activity  product   (AP)    to  the 
saturation   constant  (KT).     The  AP  of  an  undersaturated  water  is   less 
than  one,  so  log   (AP/KT)   is  negative.     Conversely,   log   (AP/KT)   of  an 
oversaturated  water  is   positive. 
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In   the  previous   general    discussion  of  carbonate  chemistry, 
equations   2,   3,   and  4  were  given  showing  the  interreactions   among  the 
carbonate  species  CO  ^  (gas),   HCO3",   and  CO3"  ^     According  to  these 
equations,   a  change  in  any  one  of  the  carbonate  species  will    produce  a 
change  in   the  others.     Combining  equations   2,   3,  and  4,  for  example, 
gi  ves 

2HC0  3"->C0  2(gas)  +  C0  3-2+H  ^0, 

which  shows    that  i f  CO  2  gas   leaves  a  water,   responding  to  a  decreasing 

P  -  2 

CO  2,   the  CO3        concentration  may  increase.     This   reaction  in  turn  will 

increase   the  AP  for  calcite  and  thus   the  degree  of  saturation  of  the 

water  with   respect  to  calcite. 

Certain  of  the  hot  springs   visibly  evolve  gas  where  they  emerge  at 

the  earth's   surface,  and  because   this   gas  is   in  part  CO  2  ,   the  carbonate 

chemistry  of  individual   springs   should  differ  as   a   function   of  the  amount 

of  gas  evolved.     Measured  and  analyzed  CO^  partial   pressures  (  ^^CO^j  from 

some  of  the  springs   are  shown  in  figure  9.     Comparison  of  this   figure 

with  the  total   dissol  ved-carbonate  data  in  table  9  shows   that  a  spring 

such   as   number  23  that  has   a  low  ^C02  value,  which  presumably   results 

from  evolution   of  a  significant  amount  of  dissolved  gas,  has   a  lower 

total    dissol ved-carbonate  concentration  than  such  springs   as  numbers  42 

and  49  that  have  higher  ^00  2  values.     The  ^CO     of  these  waters  is  as 

high   as   0.035  atmosphere   (fig.  9);   whereas,   the  ^CO.   of  the  earth's 

atmosphere   to  which   they   are  exposed  at  the  surface  is   only  0.0003 

atmosphere.      The  hot-springs  waters   certainly   could  lose  CO,,  to  the 

atmosphere  as   they  approach   the  surface. 
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To  follow  the  influence   of  CO2   outgassing  on  the   overall    carbonate 
chemistry  of  the  soring  waters,   calculations   of  log   (AP/KT)   at   various 
values  of  ^CO^  were  made   for  a   relatively  cocl   spring   (no.    23,   56.2°C) 
and  for  a  hot  spring   (no.    49,   61.8''C).      The   results  of  these  calculations 
are  shown   in  figure   11.      The  path   of  decreasing  calcite  saturation  with 
increasing     COn   for  spring  number  23  is  slightly  offset  from  that  of 
number  49,  because  of  the   temperature   differences   between   the  springs. 
More   important,   the  points   corresponding   to  the  other  springs   sampled 
lie   on  or  between  the   two   calculated   lines.      From  this,  we   conclude 
that   the   differences   among   the  hot  springs   in   pH  and  HCO3"    (table  6), 
in  total    dissolved  carbonate   (table  9),  and  in     CO^   (fig.   9)   are  super- 
ficial,  and  due  only  to  different  amounts  of  outgassing  in   the  near- 
surface  distributary  system  to   the  springs. 

Figure   11    shows    that  certain   of  the   springs   are  oversaturated 
with   respect  to  calcite.     Table  6,   though,  shows   that  the  calcium 
concentration  of  the   springs   is    identical    (44.5±0.5  mg/1),   and  that 
the  springs   that  were  sampled  had  not  yet  begun  to  deposit  calcite, 
though   they  have  the  potential    to  do  so.      If  allowed  to  flow  freely 
on   the  surface,   all    the  springs  would  continue  to   lose  CO2   and  thus 
would  deposit  calcite.     The  tufa  associated  with  the  springs  shows 
such   deposition   does   occur. 

One  significant  difference  between   the   hot-springs  waters   and 
the  cold  wells   and  springs   in   the  region  is  their  sodium  (Na)   concentration. 
The  average  sodium  concentration  of  the   12   relevant  cold-water  samples 
is   1.8  mg/1    (0.08  millimoles   per  liter),  whereas   that  of  the  hot  springs 
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Figure  11.  — Relation  between  calcite  saturation  and 
PCO^    in    the   hot  sprinq'. . 
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is  4.0  mq/1  (0.18  millinoles  per  liter).  There  is  no  difference  between 
the  carbonate,  sulfate,  or  chloride  concentrations  of  the  cold  and  hot 
springs  to  balance  the  additional  sodium  as  there  would  be  if  the  sodium 
were  a  product  of  the  solution  of  some  mineral.   Instead,  the  sodium  must 
result  from  a  reaction  by  which  it  is  exchanged  for  some  other  positively 
charged  dissolved  ion.  Such  a  reaction  can  be  written 

Ca'^-'+Na2X-CaX+2Na"^,  (13) 

in  which  Ca  "  represents  an  alkaline  earth  and  Na^X  and  CaX  represent 
solids,  perhaps,  clay  minerals,  in  the  aquifer  system,  with  X  the  sub- 
strate on  which  the  exchange  occurs. 

This  exchange  reaction  must  be  considered  in  calculating  the  pro- 
Dortion  of  plant  to  total  carbonate  to  adjust  the  measured-^'^C  concentra- 
tions of  the  hot  springs  and  to  determine  their  water  ages.  To  calculate 
this  proportion  for  the  cold  springs,  the  expression 

P=l_(  ^alkaline  earth  \  (10) 

\  total  carbonate/ 
was  used.  In  deriving  equation  10,  it  was  assumed  that  only  mineral 


carbonate  solution  affected  the  alkaline-earth  concentration  of  the  waters, 
and  the  validity  of  this  assumption  was  borne  out  by  the  agreement  between 
the  calculated  and  measured  stable  carbon-isotope  concentration  (6^^C)  of 
the  samples   In  the  hot-springs  waters,  the  alkaline-earth  concentration 
is  dependent  on  exchange  reaction  13,  as  well  as  on  the  carbonate 
chemistry,  thus  equation  10  cannot  be  used. 
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In  discussing  the  carbonate  chemistry  of  the  cold  springs,  the 
relationship  between  the  stable  carbon-isotope  concentration  (s^^C)  of  a 
sample  and  the  proportion  of  plant  to  mineral  carbonate 

613C     ,  =6l^C  ,  J^^^^^  (11) 

sample    plant' C^q^^^  /  v  ^ 

was  presented.  This  equation  can  be  used  directly  to  find  the  ^'re- 
adjustment factor,  P: 

V^tOtal/    ^i'Cp^^nt 

By  using  an  average  5^  3Cp-|gp|-=24%-.  as  before  (eq.  12),  this  becomes: 

''=«''C3,,„^.,3/-24^:  .  (14) 

Isotopic  fractionation  takes  place  between  C0„~^,  HCO^") and  CO^lgas) 

-2  _ 

as  outgassinq  and  CO3   production  from  HCO   occur.  However,  at  the 
temperature  of  these  springs  (^62'X)  the  fractionation  accompanying 
reaction  HCO3  ^H  ->C02+H20  is  nearly  equal  to  and  opposite  that  of  the 
reaction  HC03~^-H  +003"''  (Pearson,  Bedinger,  and  Jones,  1972).  Thus, 
isotopic  fractionation  will  not  interfere  with  the  calculation  of  P  by 
equation  14. 

Because  the  hot  springs  are  chemically  and  isotopically  so  similar, 
their  average  adjusted-i'+C  concentrations  and  ages  were  calculated 
(table  9).  Only  those  samples  for  which  both  ^^C  and  I'+C  data  are 
available  were  used  in  this  average.  As  table  9  indicates,  the  proportion 
of  plant  carbonate  to  total  carbonate  in  them,  from  equation  14,  is  0.608, 
and  their  adjusted-i'+C  concentration  is  58.5-percent  modern,  corresponding 
to  a  water  age  of  4,430  years. 
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Although  the  total  carbonate  concentration  of  the  hot-springs  water 
before  near-surface  nas  evolution  is  not  known,  it  certainly  cannot  be 
lower  than  the  highest  value  measured— 3.18  millimoles  per  liter.  This 
value  is  plotted  with  the  adjusted  age  on  figure  12  for  comparison  with 
the  ages  and  carbonate  concentrations  of  the  cold  wells  and  springs. 
That  the  hot  springs,  cold  sprinas,  and  wells  are   part  of  the  same  group 
is  clear,  and  is  further  evidence  that  all  are  part  of  the  same  hydro- 
logic  system. 

A  hydrologically  important  inference  can  be  drawn  from  figure  12. 
The  hot  springs  are   on  a  trend  line  defined  by  both  older  and  younger 
cold-water  samples.  Thus,  it  is  probable  that  the  age  of  the  hot-springs 
water  is  due  largely  to  its  residence  time  in  the  cold-water  part  of  the 
flow  system,  rather  than  to  an  extended  period  in  a  zone  of  heating.  A 
rapid  travel  time  through  the  heated  part  of  the  system  is  in  keeping 
with  earlier  conclusions  from  oxygen-isotope  (^^0)  exchange,  suggesting 
no  prolonged  period  of  heating,  and  the  fact  that  the  surface  temperatures 
of  the  hot  springs  closely  approach  their  probable  maximum  temperatures, 
as  calculated  from  their  silica  concentrations,  implies  rapid  flow  which 
prevents  significant  heat  loss  as  the  heated  water  rises  to  the  surface. 

Radioactivity  of  the  Waters  in  the  Hot-Springs  Region 

The  radioactivity  of  the  waters  of  the  hot  springs  has  been  studied 
by  several  investigators.  The  first  study  was  by  Haywood  (1902),  followed 
by  Boltwood  (1905),  Schlundt  (1935),  and  Kuroda  and  others  (1953).  Much 
of  the  early  interest  was  because  of  the  balneological  use  of  the  water. 
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In  1953  studies  by  Arndt  and  Damon  were  sponsored  by  the  U.S.  Atomic 
Energy  Commission,  whose  interest  was  the  radon  concentration  and  the 
source  of  the  radioactivity. 

The  presence  of  radium  in  the  waters  of  the  hot  springs  was 
established  by  Schlundt  (1935)  when  he  determined  an  average  value  of 
1.38  picograms  per  liter  (1.35  picocuries  per  liter)  of  radium  for  three 
samples.  Recent  (July  1973)  analyses  of  the  water  by  the  U.S.  Environ- 
mental Protection  Agency  show  a  radium  concentration  of  2.li0.22  pico- 
curies per  liter.  Radon,  a  gas  that  is  a  radioactive  decay  product  of 
radium,  was  analyzed  in  the  hot  springs  by  Boltwood  (1905),  Schlundt 
(1935),  and  Kuroda  and  others  (1953).  Boltwood  (1905)  reported  the  radon 
concentration  of  45  springs,  ranging  from  0.017  to  9.03  nanocuries  per 
liter  (10~^  curie/1),  with  a  modal  value  of  0.466  nanocurie  per  liter, 
Schlundt  found  the  radon  concentration  of  six  hot  springs  to  range  from 
0.125  to  0.46  nanocurie  per  liter,  with  a  modal  value  of  0.32  nanocurie 
per  liter.  Kuroda  and  others  (1953)  reported  the  radon  concentration  of  25 
hot  springs  to  range  from  0.14  to  30.5  nanocuries  per  liter,  with  a  modal 
value  of  0.82  nanocurie  per  liter. 

The  variations  from  spring  to  spring  found  by  each  investigator 
reflect  the  combined  influence  of  differences  in  analytical  methods  and 
natural  variations. 

The  source  of  radium  and  radon  in  the  hot-springs  waters  is  not 
definitely  known,  nor  has  the  water  been  analyzed  for  the  presence  of 
other  radioactive  elements.  The  presence  of  radium  and  radon  in  the 
region  is  not  peculiar  to  the  water  of  the  hot  springs.  Waters  from 
deep  wells  at  Hope,  68  miles  (109.4X10^  m)  southwest  of  Hot  Springs, 
and  at  Prescott,  53  miles  (85.3X10"'  m)  southwest  of  Hot  Springs, 
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ranged  in  radon  concentration  from  0.05  to  1.88  nanocuries  per  liter 
(Kuroda,  1953).  Waters  from  the  warm  and  cold  springs  near  Caddo  Gap, 
33  miles  (53.1X10^  m)  west  of  Hot  Springs,  ranged  in  radon  concentration 
from  0.13  to  1.85  nanocuries  per  liter  (Kuroda,  1953).  Cold  springs 
in  the  immediate  vicinity  of  Hot  Springs  also  contained  radium  and 
radon.  The  U.S.  Environmental  Protection  Agency  (written  commun., 
July  1973)  reported  a  radon  concentration  of  Happy  Hollow  Spring  of 
1.0-0.15  picocuries  per  liter  and  a  radon  concentration  of  Whittington 
Avenue  Spring  of  0.58+0.12  picocuries  per  liter.  Kuroda  (1953)  reported 
the  following  radon  concentrations,  all  in  nanocuries  per  liter: 
Whittington  Avenue  Spring,  0.36;  Whittington  Park  well,  0.03;  Happy 
Hollow  Spring,  0.74;  and  Sleepy  Valley  Spring,  4.37  and  2.80.   In 
addition,  Kuroda  (1953)  reported  the  radioactivity  of  four  samples  of 
rainwater  ranging  from  2.72  to  6.37  nanocuries  per  liter.  Radon  concentra- 
tion of  spring  waters  issuing  near  the  uranium-vanadium-niobium  bearing 
deposits  at  Potash  Sulphur  Springs,  6  miles  (9.65X10^  m)  southeast  of 
Hot  Sorings,  averaged  about  15  nanocuries  per  liter  (Arndt  and  Damon, 
1953). 

THE  HOT-SPRINGS-FLOW  SYSTEM 

The  geologic,  hydrologic,  and  water-chemistry  data  and  interpretations 
presented  in  previous  sections  of  the  report  provide  the  bases  from  which 
a  conceptual  model  of  the  hot-springs-flow  system  can  be  drawn.  Digital 
models  of  the  flow  system  provide  a  test  of  the  conceptual  model  and  a 
vehicle  for  further  elaboration  on  a  feasible  conceptual  model. 
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Conceptual  Model 

The  hot-springs  water  is  meteoric;  that  is,  it  is  derived  from 
precipitation.  The  water  is  recharged  to  formations  in  the  hot-springs 
region  (within  a  few  tens  of  miles),  as  opposed  to  being  recharged  to 
formations  distant  (several  tens  or  hundreds  of  miles)  from  the  hot 
springs.  The  origin  and  proximity  of  recharge  of  the  hot-springs  water 
are  revealed  by  the  chemical  constituents  and  isotopes  of  the  water  and 
the  flow  variations  of  the  hot  sprinas. 

The  formation  or  formations  mat  form  the  recharge  area  must  possess 
the  following  general  characteristics: 

1.  The  outcrop  area  must  be  relatively  large  and  permeable, 

that  is,  from  3  to  10  square  miles  (2.6X10-'  to 
7.8X10''  m^)  assuming  permeabilities  that  would  permit 
average  recharge  rates  of  from  2  to  6  inches 
(5.1X10"^  to  1.5X10"^  m)  per  year. 

2.  The  elevations  of  water  levels  in  the  recharge  area  must 

be  high  enough  to  provide  hydraulic  head  for  spring 
flow.  The  springs  emerge  at  elevations  ranging  from 
576  to  683  feet  (175.6  to  208.2  m)  above  mean  sea  level. 

3.  The  recharge  formations  must  be  hydraulical ly  connected  to 

the  permeable  zones  that  feed  the  springs. 
Consideration  of  the  lithology  and  structure  of  rocks  in  the  area 
reveals  two  formations  whose  outcrops  may  serve  as  recharge  areas-— the 
Bigfork  Chert  and  the  Arkansas  Novaculite.  The  outcrop  area  of  the 
Bigfork  Chert  meets  several  of  the  requirements  of  the  hot-springs 
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recharge  area — The  Bigfork  crops  out  over  an  area  of  about  36  square 
miles  north  and  northeast  of  the  hot  springs,  in  the  area  shown  in 
figure  7.  The  Bigfork  generally  possesses  moderately  high  fracture 
and  intergranular  permeability.  Faults  and  fracture  zones  are  common. 
Permeable  zones  from  the  Bigfork  Chert  to  the  hot  springs  could  be  pro- 
vided by  faults  (fig.  3)  that  juxtapose  the  Bigfork  Chert  with  permeable 
zones  in  the  Arkansas  Novaculite  and  the  Hot  Springs  Sandstone  or  provide 
a  permeable  fault  zone  through  the  Missouri  Mountain  and  the  Polk  Creek 
Shales  to  the  hot  springs. 

The  Arkansas  Novaculite  crops  out  in  a  smaller  area  in  the  vicinity 
of  the  hot  springs  than  does  the  Bigfork  Chert  (fig.  7).  Within  the 
mapped  area  in  figure  7,  the  outcrop  of  the  Arkansas  Novaculite  covers 
about  13  square  miles  (33.7X10^  m^).  The  permeability  of  the  Arkansas 
Novaculite  differs  greatly;  in  places  it  is  moderately  high  in  permeability. 
The  outcroD  of  Arkansas  Novaculite  occupies  the  highest  elevation  in  the 
vicinity  of  the  hot  springs. 

The  elevation  of  head  in  the  recharge  area  required  to  provide  flow 
depends  on  such  factors  as  transmissi vity  of  the  aquifer,  areal  extent  of 
the  recharge  area,  distance  of  the  recharge  area  from  the  springs,  depth 
of  flow  of  water  in  route  to  the  springs,  and  the  temperature  of  the  water 
in  the  aquifer.  These  factors,  except  for  the  effect  of  temperature  on 
hydraulic  head,  are  accounted  for  in  the  digital  models  of  the  aquifer. 
The  effect  of  the  increased  temperature  of  water  in  the  flow  system  reduces 
the  head  required  in  the  recharge  area  to  drive  the  system.  The  density 
of  water  above  4.0°C  (39.2"F)  decreases  as  its  temperature  increases. 
Water  at  61.7''C  (143. 0°F)  (average  hot-springs  temperature)  is  98  percent 
as  heavy  as  water  at  17.7°C  (63.9°F)  (average  temperature  of  recharge  water). 
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Therefore,  a  column  of  water  at  17.7'^C  (63.9"F),  1,000  feet  in  length, 
will  support  a  1,020-foot  column  of  water  at  61.7°C  (143. 0°F).  It  would 
be  possible  to  have  a  thermal  artesian  system  in  which  the  water-level 
elevation  in  the  recharge  area  is  less  than  the  water-level  elevation 
in  the  discharge  area. 

Water  levels  in  ^'he  outcrop  area  of  the  Bigfork  Chert  range  from  a 
few  feet  above  land  surface  (flowing  wells)  to  30  feet  (9.14  m)  below 
land  surface.  Water  levels  in  the  Arkansas  Novaculite  range  from  a  few 
feet  above  land  surface  to  50  feet  (15.2  m)  below  land  surface.  As  an 
approximation  of  the  area  in  which  heads  in  the  Bigfork  Chert  and  the 
Arkansas  Novaculite  are  sufficient  to  drive  the  artesian  system,  the 
outcroD  area  above  700  feet  (213.4  m)  mean  sea  level  was  chosen.  Within 
the  area  shown  in  figure  11,  about  8,5  square  miles  (22.0X10^  m^)  of  the 
outcrop  of  Bigfork  Chert  is  above  700  feet  (213.4  m)  mean  sea  level 
and  almost  all  of  the  outcrop  of  Arkansas  Novaculite  is  above  700  feet 
(213.4  m)  mean  sea  level.  However,  under  about  5  square  miles  (13.0X10^  m^) 
of  the  Arkansas  novaculite  outcrop,  the  northwest  limit  of  the  anticline 
north  of  the  hot  sprinqs,  the  formation  dips  northwestwardly  and 
probably  is  not  in  a  favorable  structural  position  to  supply  water  to  the 
hot  springs. 

Analysis  of  streamflow  provides  information  on  the  magnitude  of 
recharge.  Stream  discharge  and  rainfall  records  have  been  collected  on 
three  small  basins  in  the  vicinity  of  the  hot  springs.  The  locations  of 
the  gages  on  East  Fork  and  West  Fork  of  Hot  Springs  Creek  are  shown  in 
figure  7.  The  gage  on  Glazypeau  Creek  is  located  about  9  miles  (14.5X10^  m) 
north-northeast  of  the  hot  springs.  Each  basin  is  relatively  small  and 
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is  underlain  principally  by  Bigfork  Chert.  The  Arkansas  Novaculite 
underlies  the  higher  elevations  of  East  Fork  and  West  Fork  of  Hot 
Springs  Creek  (fig.  7). 

Base  flow  is  ground-water  discharge  to  the  stream.  This  part  of 
the  streamflow  is  derived  from  recharge  of  precipitation.  The  base-flow 
figures  (table  10)  show  that  a  relatively  large  part  of  rainfall  is 
recharged.  These  figures  are  primarily  a  measure  of  the  recharge  to  the 
Bigfork  Chert.  Also,  the  base  flow  indicates  that  interbasin  transfers 
of  ground  water  occur.  Such  interbasin  transfers  of  ground  water,  probably 
from  several  stream  basins,  supply  the  water  to  the  hot  springs. 

The  large  surface  area,  the  great  capacity  to  admit  recharge,  and 
the  relatively  high  permeability  of  the  Bigfork  suggest  that  the  Bigfork 
Chert  is  the  principal  recharge  source  to  the  hot-springs  artesian  system. 
The  Arkansas  Novaculite  has  a  smaller  outcrop  area  of  potential  recharge, 
a  lower  capacity  to  admit  recharge,  and  a  generally  lower  permeability. 
The  Arkansas  Novaculite  is  probably  a  source  of  part  of  the  recharge  to 
the  hot-springs  system.  Further  analysis  of  the  Bigfork  Chert  and 
Arkansas  Novaculite  as  the  recharge  sources  and  principal  aquifers  of 
the  flow  system  is  made  in  the  section  on  modeling. 

The  ^^C  age  of  the  hot-springs  water  averages  4,430  years.  The 
greater  part  of  the  time  that  the  water  is  in  underground  circulation 
the  movement  is  very  slow.  That  is,  movement  is  in  the  order  of  a  few 
feet  to  a  few  tens  of  feet  per  year  and  is  distributed  through  a  large 
volume  of  the  aquifer.  This  slow  movement  of  water  continues  to  great 
depth  in  the  aquifer,  where  the  water  is  heated  by  contact  with  rocks  of 
high  temperature.  The  depth  of  circulation  is  not  known,  but  could  be 
from  2,000  to  8,000  feet  (610  to  2,400  m). 
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The  high  temperature  of  the  rocks  at  depth  has  been  attributed 
to  (1)  the  normal  geothermal  gradient,  which  is  about  0.01°C  (O.OIS^F) 
per  foot  (0.3048  m),  and  (2)  an  abnormally  high  geothermal  gradient. 
An  abnormally  high  thermal  gradient  could  be  created  by  a  radioactive 
or  an  igneous  mass. 

Highly  permeable  zones,  probably  related  to  jointing  or  thrust 
faulting,  collect  the  heated  water  in  the  aquifer  and  provide  avenues 
for  the  water  to  travel  to  the  surface.  Rapid  movement  of  the  water 
from  depth  to  the  surface  is  indicated  by  the  very  small  decrease  in 
temperature  from  the  maximum  temperature  attained  at  depth  to  the  surface. 

Digital  Flow  Models 

The  purpose  of  modeling  the  hot-springs-flow  system  was  to  test 
several  hypotheses  of  the  flow  system.  Each  hypothesis  requires  details 
of  aquifer  geometry  and  various  combinations  of  parameters,  such  as  area 
and  depth  of  circulation,  hydraulic  conductivity,  and  porosity.  The 
model  analysis  provides  solutions  to  parameters  that  can  be  compared  with 
observed  data  for  the  hot  springs  (such  as  temperature,  discharge,  ^^C 
concentration,  and  silica  concentration).  Two  physically  different  models 
were  applied  in  the  analysis  of  the  hot  springs,  one  a  steady-state 
hydrothermal  model  consisting  of  three  parts— a  heat-flow  model,  a  ground- 
water-flow  model,  and  a  group  of  routines  to  calculate  variables  (silica, 
time  of  travel,  ^''C)  that  are   functions  of  the  ground-water-flow  and 
heat-flow  patterns.  This  steady-state  hydrothermal  model  is  described 
in  the  following  paragraphs.  The  other  model,  a  nonsteady  ground-water- 
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flow  model,  was  used  to  simulate  the  transient  response  of  spring  flow 
to  seasonal  variations  in  recharge.  The  nonsteady  ground-water  model 
is  described  by  Bedinger,  Reed,  and  Griffin  (1973). 

Description  of  the  Hydrothermal  Model 

The  heat-flow  and  ground-water-flow  models  are  based  on  the  equations 
of  continuity  for  heat  and  water,  respectively.  That  is,  the  net  flux 
of  heat,  or  water,  into  an  element  of  volume  is  equal  to  the  change  in 
heat,  or  water,  concentration  of  the  element.  Both  these  models  are  steady- 
state  models;  it  is  assumed  that  neither  temperature  nor  head  change  with 
time.  Therefore,  the  net  flux  of  heat,  or  water,  into  any  element  of 
the  model  is  equal  to  zero.  The  equation  of  continuity  for  heat,  or 
water,  is  written  for  each  of  the  elements  of  the  model.  If  there  are 
N  elements  in  the  model,  the  result  is  N  linear  equations  in  I^  unknown 
temperatures,  or  heads.  These  N  linear  equations  are   solved  by  (1)  assuming 
an  initial  temperature,  or  head;  (2)  using  the  Jacobi  iteration  method 
to  calculate  new  values  for  temperature,  or  head;  and  (3)  calculating 
the  residual  for  each  element.  The  residual  is  the  sum  of  the  head,  or 
water,  inflow  to  each  element.  The  residual  would  be  zero  for  each 
element  if  the  N[  temperatures,  or  heads,  were  exact  solutions  of  the  N_ 
equations.  The  initial  values  used  were  18.3°C  (65.0°F)  for  the  heat- 
flow  model  and  0  feet  head  for  the  ground-water-flow  model.  Steps  2  and 
3  are  rppeated  until  the  sum  of  all  tho  residuals  in  the  model  is  less 
than  0.1  percent  of  the  total  flux  in  the  model.  The  total  flux  is  equal 
to  the  total  heat  inflow  to  the  heat-flow  model,  or  to  the  total  recharge 
for  the  ground-water-flow  model. 

76 


The  equation  for  the  heat-flow  model  is: 

H=.|((x,.  )(-.;^>(vx„)(M.)).  IE((VXJ  (!^ 

+(x.-Xo)('^'^-^    SP  ^iQi+X,0.-HX3Q3+X,Q,-XoQo)  . 

The  subscripts  in  the  above  equation  refer  to  node  location  (a  node  is 
the  center  of  a  volume  element);  o  refers  to  the  nodal  location  being 
considered,  i  refers  to  the  node  above,  :.  refers  to  the  node  to  the 
right,  3  refers  to  the  node  below,  and  4  refers  to  the  node  to  the  left. 
In  the  above  equation,  H  is  the  heat-flow  residual,  in  BTU/day;  KN  is 
the  thermal  conductivity,  in  BTU/day. ft. °F;  X  is  the  temperature,  in 
^F;  XD  is  the  horizontal  node  spacina,  in  feet;  YD  is  the  vertical  node 
spacing,  in  feet;  MA  is  the  thickness  of  the  aquifer  (equivalent  to 
thickness  of  model),  in  feet;  SP  is  the  volumetric  specific  heat  of  water, 
62.4  BTU/°F-ft3  (6.7X10'^  joule/(°C.m3) ;  Q,  ,  .  .  is  the  inflow  of  water 
from  node  1,  2,  3,  and  4,  respectively,  in  ft^/day;  and  Q  is  the  outflow 
of  water  from  the  element,  in  ftVday. 

The  equation  for  the  nround-v/ater-flow  model  is: 

QT  is  the  ground-water-flow  residual,  in  ft^/day;  P  is  the  hydraulic  con- 
ductivity, in  ft/day;  and  S  is  the  ground-water  potential  (head),  in  feet. 
Other  symbols  and  the  subscripts  are  as  explained  for  the  heat-flow  model. 
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According  to  R.  0.  Fournier  (oral  commun.,  1972),  the  relation 
between  saturated  silica  concentration  and  temperature  is 


MGL  =  60,060  (10-;  ^^—^  +   0.09/'. 

where  MGL  is  the  silica  concentration  at  saturation  with  chalcedony,  in 
mg/1;  and  K  is  the  temperature,  in  degrees  Kelvin  (=°C  +  273.2).  The 
assumptions  made  in  the  model  of  silica  concentration  are:  (1)  If  the 
water  is  supersaturated  with  silica  (silica  concentration  greater  than 
that  given  by  the  above  equation),  the  silica  will  remain  in  solution 
and  not  precipitate;  (2)  if  the  water  is  undersaturated  with  silica 
(silica  concentration  less  than  that  indicated  by  the  above  equation), 
sufficient  silica  will  go  into  solution  so  that  the  water  becomes  saturated 
with  silica.  The  effect  of  these  two  assumptions  requires  that  the  water 
either  be  saturated  or  supersaturated  with  silica.  The  iterative  procedure 
for  calculating  the  silica  distribution  in  the  aquifer  is  as  follows. 
First,  the  silica  concentration  at  a  node  was  calculated  from  the  inflow 
by 


SI 


QjSIL^+Q,SIL.+Q,SIL,iQ,SIL, 


Q  +Q  +Q  -HQ 

^1  V.   3  ^if 


where  SIL  is  the  silica  concentration  at  a  node  as  calculated  from  the 
inflow,  in  mg/1;  SIL       are  silica  concentrations  at  surrounding 
nodes,  in  mg/1;  and  Q  ^    are  inflows  (zero,  if  outflow)  from  sur- 

1  ,  ^  ,  3  ,  H 

rounding  nodes,  in  ft-^/day.  Subscripts  are  as  defined  in  the  heat-flow 
equation.  The  silica  concentration  at  the  node  is  chosen  as  the  larger 
of  SIL,  given  above,  or  MGL,  as  given  in  the  previous  equation.  On  the 
first  iteration,  all  the  silica  values  are  as  given  by  the  equation  for  MGL, 
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Succeeding  iterations  produce  silica  values  as  weighted  from  inflow  by 

the  equation  for  SIL.  An  indicator  is  kept  for  the  naximum  change  in 

silica  fron  one  iteration  to  the  next.  When  this  indicator  becomes 

less  than  0.0001,  the  final  silica  values  have  been  computed  and  the 

iteration  is  stopped. 

The  iterative  procedure  for  calculating  time  of  travel  uses  the 

equation 

.,^  ...  jX!l)jYnHPORi 

Qi+q7q".-^+  (B)(xn){vul       ' 

(Dj(MAj+MA,,)-;D   (MA  +MA  J-t-n   ("A  +KAj,)+DjMA^+!-1Ag))    , 

where  TF  is  the  ti-^^  of  travel,  in  years;  D  is  one  if  the  corresponding 

Q  is  positive  and  zero  if  the  corresDonding  Q  is  zero;  POR  is  the  porosity, 

dinensionless ;  and  b  is  the  recharge  rate  at  the  node,  in  ft/day.  Other 

symbols  and  the  subscripts  are  as  defined  previously.  An  indicator  is 

used  that  is  the  maximum  change  in  time  of  travel  from  one  iteration  to 

the  next.  When  this  indicator  becomes  less  than,  or  equal  to,  0.0001, 

the  iteration  is  stopped. 

The  model  for  ^^C  is  based  upon  the  exponential  decay  of  ^^C,  a 

ratioactive  isotope  of  carbon.  This  decay  can  be  expressed  as 

CU:t=C14e-^^-'057XI0-")t, 

where  C14  is  t^""  initial  ^'^"^   ^^ntent,  C14^  is  the  ^^C  concentration  after 
a  period  of  time  has  elapsed,  and  t  is  the  elapsed  time  in  years.  The 
recharge  area  of  the  aquifer  is  not  the  only  source  for  the  carbon  in 
the  water,  as  some  carbon  is  added  by  solution  of  aquifer  material.  As 
table  9  shows,  within  the  bulk  of  the  system,  only  55  to  65  percent  of 
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the  carbonate  present  is  ^^C  bearing.  The  ^^^C  model  uses  65  percent 
as  an  initial  value  and  decreases  this  value  accordinq  to  the  exponential 
decay.  The  ^'^C  model  is  an  iterative  procedure  calculating  ^"^C  from 
the  equations 


r 


14 


_QiC14^e  ^fQ.CU.e  '+Q3C143e^^+Q„C14,eSB(XD)(YD) 
~^  Qj+Q2+Q3+Q^+B(XD)(YD) 

and  ^ 

(-1.2097X10-  )(XD)(YD)(DOR)(MA  +MAo) 

'^-   ' mT, ' '  ^"^'^'^'^  ' 

where  Cl-^  is  the  ^^C  concentration,  in  percent  of  modern  as  measured, 
and  other  symbols  and  subscript  meanings  are   as  used  previously.  An 
indicator  is  used  whose  value  is  the  maximum  change  of  CI 4  from  one 
iteration  to  the  next.  When  this  indicator  becomes  less  or  equal  to 
0.0001,  the  iteration  is  terminated. 


Aoplication  of  Dinital  Models  to  the  Hot-Sprinqs-Flow  System 

Most  of  the  data  on  the  hot-springs-flow  system,  such  as  water 
temperature,  water  samples  for  analysis,  and  flow  measurements,  have 
been  collected  from  points  at  or  near  the  earth's  surface.  Subsurface 
data  are  from  wells,  the  deepest  of  which  is  about  335  feet  (102.4  m). 
Conditions  at  depth  must  be  inferred  by  such  means  as  projection  of 
surface  dips  of  rock  formations  and  faults,  and  changes  in  temperature 
and  chemical  nature  of  the  water  from  the  recharge  area  to  the  springs. 
The  digital  models  provide  a  further  tool  with  which  conditions  at 
depth  can  be  interred.  Measured  properties  of  the  flow  system  and  the 
water  can  be  entered  into  the  digital  model.  The  flow  model  provides 

8ii 


simultaneous  solutions  to  the  equations  for  each  of  the  several  proper- 
ties  The  design  of  the  system  at  depth  can  be  changed  experimentally 
within  the  bounds  of  known  parameters  so  that  the  model  results  are 
compatible  with  all  measured  properties  related  to  the  system. 

Following  is  a  summary  of  criteria  that  must  be  specified  in 
designing  a  model  of  the  flow  system. 

1.  Thickness,  porosity,  hydraulic  conductivity,  storage 

coefficient,  and  deoth  of  the  aquifer. 

2.  Thermal  conductivity  of  the  rocks  and  temperature 

gradient  into  the  aquifer. 

3.  Recharge  to  the  aquifer. 

4.  Temperature  of  water  in  the  recharne  area. 

The  initial  flow  models  were  designed  from  projections  of  surface 
dips  of  the  aeologic  formations  and  from  known  and  estimated  hydrologic 
and  thermal  oronerties  of  the  aquifer.  Two  flow  models  were  designed 
and  tested  during  the  study.   In  one  model  the  Binfork  Chert  served  as 
the  principal  aquifer.  In  the  other  model  the  Arkansas  Novaculite  served 
as  the  principal  aquifer.  During  development  of  each  model,  the  design 
was  changed  in  attempting  to  achieve  harmony  between  the  input  data  and 
the  known  constraints  on  the  system. 

The  Bigfork  Chert  was  modeled  as  a  two-dimensional  olane,  inclined 
in  the  subsurface  to  the  southeast  (fig.  13).  The  flow  of  water  was 
modeled  to  a  point  in  the  subsurface  beneath  the  spring  discharge  area. 
Flow  upward  to  the  surface  was  not  modeled  because  the  model  is  limited 
to  a  two-dimensional  representation  of  the  flow  system.  Temperature  and 
silica  data  on  the  hot  springs  indicate  that  the  flow  of  water  from  depth 
to  the  surface  is  relatively  rapid.  That  the  flow  from  depth  to  the 
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Figure  13.— Model  of  flow  in  the  Bigfork  Chert. 
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surface  was  not  modeled  is  not  considered  significant  with  regard  to 

conclusions  drawn  from  the  model.  A  summary  of  the  input  data  for  this 

model  is  shown  in  the  following  table: 

Thickness  of  aquifer feet--        1,500 

Porosity 0.20 

Hydraulic  conductivity  at  18.3°C 

(nonuniform) ft/day--       0.5-10 

Thermal  conductivity cal/sec-cm-°C--       0.0062 

Temperature  gradient  (nonuniform) °c/ft--    0.011-0.11 

Recharge  rate  (nonuniform) ft/yr--     0.05-0.25 

The  Bigfork  Chert  in  the  model  is  thicker  than  the  average  recon- 
structed stratigraphic  thickness  for  this  formation.  However,  the  Big- 
fork  Chert  contains  multiple  folding  which  increases  the  actual  thickness 
of  the  water-bearing  zone.  Also,  the  underlying  Womble  Shale  contains 
permeable  zones  of  limestone  and  chert  which  effectively  increase  the 
thickness  of  the  contiguous  water-bearing  zone.  The  hydraulic  conductivity 
of  the  Bigfork  Chert  was  modeled  as  1  ft/day  (3.53X10"^  m/s)  over  the  out- 
crop area  and  in  most  of  the  subsurface.  The  higher  values  of  hydraulic 
conductivity  were  modeled  in  a  zone  trending  east-northeast  from  the  hot 
springs,  parallel  to  the  trace  of  the  thrust  faults  through  the  hot-springs 
area  (fig.  2).  The  hydraulic  conductivity  in  the  model,  after  being 
adjusted  by  the  program  to  reflect  the  temperature  of  the  water,  ranged 
from  1.  to  25.  ft/day  (3.53X10"^  to  8.82X10"^  m/s).  The  maximum  depth  of 
water  circulation  in  the  model  is  6,000  feet  (1,830  m).  With  this  depth 
of  circulation  and  a  thermal  conductivity  of  0.0062  cal/sec-cm-°C,  a  non- 
uniform temperature  gradient  of  from  normal  (0.011°C/ft  (3.61X10"^  °C/s)) 
to  O.irc/ft  (3.61X10'i  °C/s)  was  specified  in  the  model  in  order  to 
achieve  the  observed  temperature  at  the  springs.  Results  of  the  model 

are  summarized  in  table  10. 
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Table   10. — Comparison  of  the  Bigfork  Chert  modeled  and  observed 

data  for  the  hot  springs 

Parameter  Modeled    Observed 


Temperature °C--  61.2  61.7 

Flow ftVday—  0.97X105  1.1X10^ 

Silica mg/1—  40.7  41.7 

Carbon-14  age years—  5,026  4,430 

Maximum  head  in  recharge  area 

above  hot  springs feet —  18                    260 
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The  results  of  the  Bigfork  Chert  model  indicate  that  the  known  and 
assumed  parameters  fit  well  within  the  constraints  on  the  system.  The 
model  indicates  that  the  available  head  difference  is  ample  to  provide 
the  hydraulic  pressure  for  the  spring  flow.  An  additional  component 
of  head  is  produced  by  the  decreased  density  of  the  water  at  tempera- 
tures higher  than  the  assumed  temperature  in  the  recharge  area  (18.3°C 
(65.0°F)).  The  ample  head  available  to  drive  the  flow  system  provides 
a  leeway  in  several  of  the  estimated  parameters  used  in  the  model,  such 
as  thickness  and  permeability  of  the  aquifer  and  rate  and  areal  extent 
of  the  recharge.  For  example,  a  valid  model  could  be  designed  using 
lower  values  for  permeability  and  thickness  of  the  Bigfork  Chert,  or 
a  valid  model  could  be  designed  using  a  lower  recharge  rate  over  a 
larger  area  or  using  a  higher  recharge  rate  over  a  smaller  area. 

The  Arkansas  Novaculite  was  modeled  as  a  vertical  two-dimensional 
plane  (fig.  14).  Flow  was  modeled  to  a  point  in  the  subsurface  near 
the  springs  from  which  flow  to  the  surface  is  probably  rapid  and  through 
a  very  permeable  zone.  A  summary  of  input  data  for  a  model  of  the 
Arkansas  Novaculite  is  shown  in  the  following  table: 

Thickness  of  aquifer feet--        500 

Porosi  ty 0. 1 

Hydraulic  conductivity ft/day--         1. 

Thermal  conductivity cal/sec-cm-°C--      0.0062 

Temperature  gradient °C/ft--       0.011 

Recharge  rate ft/yr--       0.63 

Results  of  the  Arkansas  Novaculite  model  are  summarized  in  table  11, 
The  critical  constraint  on  the  model  of  the  Arkansas  Novaculite  is  the 
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Figure  14.— Model  of  flow  in  the  Arkansas  Novaculite. 
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Table   11. —       ■^ain.son  of  :hc  ArKansan  Uovaaulite  •^^odeled  and 
observed  en 'a  for  the  lut  svviyicr. 


Parameter  Modeled  Observed 


Temoerature C--  54.4  61.7 

Flow ft/day--  O.SBXlu'--  1.1X10^ 

Silica mq/1--  36.0  41.7 

Carbon-14  age--- years--  2,760  4,430 

Maximum  head  in   recharge  area 

above  hot  springs feet--  1,222  420 


difference  in  head  between  the  recharge  area  and  the  springs.  The  high 
head  (1,222  feet  (372.5m))  required  by  the  model  exceeds  the  available 
head  by  a  factor  of  3.  This  head  difference  is  much  larger  than  the 
head  provided  by  the  temperature  differential  (estimated  to  be  240  feet 
(73.2  m)).  Known  values  of  permeability  of  the  Arkansas  Novaculite 
indicate  that  the  permeability  cannot  reasonably  be  expected  to  be 
greater  than  1  ft/day  (3.53X10"^  m/s)  over  an  area  as  large  as  that 
modeled.  The  temperature  and  silica  comparisons  between  the  model  and 
observed  data  are  not  significant  with  regard  to  validity  of  the 
Arkansas  Novaculite  model. 

The  results  of  the  model  of  the  Arkansas  Novaculite  indicate  that 
the  Arkansas  Novaculite  is  not  likely  the  sole  aquifer  of  the  hot-springs 
system. 

A  model  of  nonsteady  flow  was  used  to  simulate  the  time  response 
of  the  spring  flow  to  seasonal  variations  in  recharge.  The  hot  springs 
have  an  annual  fluctuation  in  discharge  equal  to  about  15  to  20  percent 
of  the  maximum  flow.  The  highest  flow  occurs  in  the  spring  and  the 
lowest  flow  occurs  in  the  late  summer.  The  permeability  distribution 
used  in  the  model  was  taken  from  the  hydrothermal  model  and  included  the 
effect  of  water  temperature.  A  storage  coefficient  of  0.005  was  used 
in  the  model.   In  simulating  recharge  to  the  aquifer,  yearly  recharge 
'./as  divided  into  four  90-day  periods.  Recharge  was  distributed  areally 
in  proportion  to  the  rate  used  in  the  hydrothermal  model.  Recharge  was 
distributed  in  the  four  time  periods  as  a  factor  multiplied  by  the 
average  rate.  The  sequence  of  factors  was  1,  0.5,  1,  and  1.5. 
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MAN'S  EFFECTS  ON  THE  HOT  SPRINGS 

Man's  principal  effects  on  the  hot  springs  can  be  placed  in  three 
categories:  (1)  activities  that  reduce  the  natural  flow  of  water  in 
the  system,  (2)  activities  that  introduce  contaminants  into  the  system, 
and  (3)  activities  that  disturb  the  flow  paths  of  hot  water  in  the 
discharge  area. 

Reduction  in  Flow 

Activities  that  reduce  the  natural  flow  of  water  in  the  system 
include  the  pumping  of  ground  water  in  the  recharge  and  discharge  areas 
and  land-use  practices  that  reduce  recharge  to  the  aquifer.  A  reduction 
in  flow  in  the  system  induces  a  rapid  effect  in  spring  flow.  The 
response  of  the  system  to  changes  in  flow  is  evidenced  by  the  seasonal 
response  in  flow  to  variations  in  recharge. 

Land  Use 

Man,  by  either  design  or  accident,  exerts  considerable  influence 
on  the  amount  of  rainfall  that  enters  the  flow  system  by  such  activities 
as  lumbering;  agriculture;  quarrying;  terracing;  building  houses,  roads, 
businesses,  parking  lots,  parks,  storm  drains,  and  pumping  wells.  The 
same  activity  can  either  increase  or  decrease  recharge,  depending  upon 
the  interaction  with  the  natural  hydrologic  system. 

The  effects  of  man's  activities  to  date  on  the  hot  springs  are  not 
known  because  the  hydrologic  history  of  the  springs  — particularly  the 
history  of  the  flow  of  the  hot  springs  — is  not  known.  The  greatest 
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potential  effect  man  has  produced  on  the  spring  system  is  the  urbani- 
zation of  parts  of  Hot  Springs  Creek  and  Gulpha  Creek  basins.  Land 
use  in  these  basins,  which  comprise  the  major  part  of  the  recharge  area, 
is  shown  in  figure  15. 

Urban  land  use 


Twenty-six  percent  of  Hot  Springs  Creek  basin  and  5  percent  of 
Gulpha  Creek  basin  is  urbanized.  Urban  land  cover  in  the  basins  is 
shown  in  the  following  table. 


Cover  by   roofs,   roads,  parking 

■  ■    ■  1 

Hot  Springs 
Creek  basin 

Gulpha  Creek 
basin 

lots,   and  other  imoermeable 
material 

Urban 

area 

(sq  mi) 

Percent 
of 

urban 
area 

Urban 
area 

(sq   mi) 

Percent 
of 

urban 
area 

i^ore  than  80-percent  covered 

40-   to  80-percent  covered 

10-   to  40-percent  covered 

Less    than   10-percent  covered 

0.112 
.050 
.598 
.183 

12 

5 

64 

19 

0.0 
.014 
.165 
.015 

0 
7 

85 
8 

Total  : 

0.943 
26 

100 

0.194 
5 

100 

otjuarti    Mil  1 1;^ ------------ - 

rcrLcML     Ul      UablFI —  -- 

Analysis  of  streamflow  provides  information  on  the  effect  of  urbani- 
zation on  stream  runoff.  Stream  discharge  and  rainfall  records  have 
been  collected  on  three  small  basins  in  the  vicinity  of  the  hot  springs. 
West  Fork  of  Hot  Springs  Creek  and  Glazypeau  Creek  basins  are  principally 
rural  areas.  A  small  part  of  West  Fork  of  Hot  Springs  Creek  basin  is 
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slightly  urbanized  (fig.  15);  the  largest  part  is  forested.  Glazypeau 
Creek  is  all  rural,  covered  by  forest  land  and  pastures.  East  Fork 
of  Hot  Springs  Creek  basin  is  mostly  urbanized  (fig.  15). 

A  summary  of  rainfall  and  stream  discharge,  from  October  1970 
through  September  1971,  for  the  three  basins  is  given  in  table  12. 
Stream  discharge  was  divided  into  direct  surface  runoff  and  base  flow 
by  means  of  hydrograph  separation.  As  the  name  implies,  direct  surface 
runoff  is  water  from  rainfall  that  runs  off  over  the  land  surface  to 
the  stream  channel.  Base  flow  is  ground-water  discharge  from  the  aquifer 
to  the  stream.  Base  flow  appears  as  springs  and  seeps  along  and  adjacent 
to  the  stream  channel.  This  part  of  the  streamflow  is  derived  from 
precipitation  that  soaks  into  the  ground  and  moves  downward  to  the 
aquifer.  Thus,  base  flow  is  a  measure  of  recharge  to  the  aquifer. 

Basin  loss  is  the  difference  between  rainfall  and  total  discharge. 
Basin  loss  includes  evapotranspiration  and  gain  or  loss  by  interbasin 
movement  of  ground  water. 

The  differences  between  percentages  of  direct  runoff,  base  flow, 
and  basin  loss  of  the  stream  basins  (fig.  15)  are  functions  of  land 
cover,  topography,  other  basin  characteristics,  and  interbasin  movement 
of  ground  water.  The  reduction  in  recharge  cannot  be  measured;  how- 
ever, the  complementary  effect— that  is,  the  increase  in  direct  runoff 
of  rainfall— can  be  observed  by  comparing  the  streamflow  data  for  East 
and  West  Forks  of  Hot  Springs  Creek.  However,  any  reduction  in  recharge 
has  been  less  than  the  proportion  of  urbanization  in  the  basins.  The 
more  heavily  covered  areas  are  at  lower  elevations  and  were  probably 
not  recharge  areas  prior  to  urbanization.  Because  the  land  surface  is 
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hiHy  and  not  naturally  suited  to  dense  urbanization,  the  impermeable 
cover  in  most  of  the  urban  area  is  sparse  to  light  (less  than  40- 
percent  covered).  Thus,  the  areas  at  higher  elevations,  v/hich  are  the 
more  likely  recharge  areas,  are  less  affected  by  urbanization. 

Landscape  terracing  for  oarks,  ballfields,  trailer  parks,  lawns, 
and  other  enterprises  retards  runoff;  and  if  the  permeability  of  the 
ground  is  not  reduced  or  covered  by  an  impervious  surface,  such  ter- 
racing can  facilitate  recharge. 

Lawn  watering  might  seem  to  offer  additional  opportunity  for 
recharge.  But,  after  one  considers  that  lawns  are  generally  watered 
during  the  season  of  peak  evapotranspiration  and,  by  usual  practice, 
in  amounts  less  than  actually  required,  this  factor  is  probably 
negl igible. 

Rural  land  use 

Eighty  percent  of  Hot  Springs  Creek  and  Gulpha  Creek  basins  is 
forested  with  pine  and  deciduous  hardwood.   Forest  cover  affects  water 
yield— runoff  and  recharge— in  a  basin.  Experimental  work  has  shown 
that  water  yield  can  be  increased  by  deforestation  (Hoover,  1944). 
Timber  operations  in  the  recharge  area  can,  and  probably  do,  affect 
recharge  and  overland  runoff.  Even  timber-management  practices,  which 
can  change  the  species  makeun  of  the  forest  or  the  ratio  of  deciduous 
to  evergreen  cover,  can  affect  the  water  yield  because  of  different 
water  requirements  for  different  species  (Minckler,  1939;  and  Horton, 
1923)  and  interception  of  rainfall  by  pine  trees  is  greater  than  by 
deciduous  trees  (Davis,  1939). 
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About  4  percent  of  the  rural  area  is  deforested.  The  deforested 

area  includes  Belvedere  Country  Club,  pastureland,  and  so-called  old 
fields  which  are  reverting  to  forests.  A  very   small  part  of  the  area 

is  cultivated.  The  deforested  area  probably  permits  greater  runoff 
than  the  forested  area,  because  of  less  interception  of  rainfall  by 
foliage  and  less  transpiration  by  plants.  For  the  same  reasons,  re- 
charge may  be  greater  in  deforested  areas.  However,  recharge  in 
deforested  land  is  dependent  upon  subsequent  land  use.  Land-use 
practices  that  compact  the  soil  and  reduce  its  permeability  will  reduce 
recharge. 

Quarry  operations  orovide  accelerated  recharge  to  the  aquifer; 
but,  because  of  the  small  area  coverW  the  net  effect  is  small. 

Ground-Water  Pumping 

Pumping  wells  in  both  the  recharge  and  discharge  areas  can  have 
a  definite  effect  on  the  hot-springs  flow.  A  well  that  is  pumped  will 
draw  water  initially  from  storage  in  the  aquifer.  With  time,  the 
effect  of  pumping  will  reach  areas  of  natural  recharge  or  discharge. 
Pumping  from  wells  in  the  recharge  area  may  decrease  discharge  of  the 
cold  springs  in  the  recharge  area,  as  well  as  decrease  flow  of  the  hot 
springs  in  the  discharge  area.  There  are  few  wells  in  and  near  the 
recharge  area,  and  total  pumpage  is  probably  small.  Most  of  the  wells 
are  of  small  capacity  and  are  used  for  rural  domestic  supplies.  The 
larger  wells  include  one  at  the  Belvedere  Country  Club  and  a  well  at 
the  Majestic  Hotel . 

Pumping  from  wells  in  the  hot-springs-discharge  area  will  have 
the  immediate  effect  of  reducing  the  natural  spring  flow.  Wells  in  the 
discharge  area  are  the  emergency  water  well  at  the  Rehabilitation  Center, 
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sprina  39— a  hot  spring  enlarged  and  deepened  into  a  well  — that  formerly 
supplied  hot  water  to  the  Rehabilitation  Center,  a  well  at  the  Arlington 
Hotel,  and  two  wells  in  the  National  Park  used  for  cooling. 

Pollution 

Biological  pollution  potential  of  the  hot  springs  from  the  recharge 
area  is  not  a  problem  because  of  the  long  traveltime  and  the  great 
distance  from  the  recharge  area  to  the  springs.  Chemical  pollutants, 
on  the  contrary,  are  not  filtered  out,  although  some  chemical  changes  may 
take  place  by  reactions  between  the  chemical  pollutants  and  the  natural 
ground  water  and  the  aquifer.  Generally,  a  chemical  pollutant  entering 
the  flow  system  in  the  recharge  area  will  emerge  from  a  spring.  Threat 
of  pollution  from  chemical  contamination  by  accidental  spillage  is  re- 
mote, because  of  the  relatively  small  volume  of  pollutants  presumably 
involved  and  the  mixing  with  and  dilution  by  natural  ground  waters 
during  the  long  period  of  flow  from  the  recharge  area  to  the  springs. 
Chemical  pollution  would  be  possible  by  dumping  large  volumes  of  pol- 
lutants on  the  ground  or  in  pits  in  the  recharge  area  during  a  sustained 
period  of  time  or  extensive  and  sustained  application  of  herbicides, 
pesticides,  and  fertilizer.  Presently,  there  are   no  known  potential 
sources  of  chemical  pollution  in  the  recharge  area. 

The  hot-springs-collection  system  is  not  closed  to  surface  sources 
of  Dollution  in  the  discharge  area.  Travel  of  a  pollutant  from  the  hill- 
side above  the  spring  outlets  may  reach  the  springs  very  rapidly. 
Tritium  analyses  show  that  post-1950  water  is  a  part  of  the  spring  flow. 
Also,  the  temperature  response  of  the  Fordyce  Spring  to  heavy  precipitation 
indicates  rapid  travel  of  water  from  the  hillside  to  the  spring. 
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Flow  Disturbance  in  the  Discharge  Area 

Activities  that  disturb  the  flow  paths  of  the  hot  water  in  the 
discharge  area  include  excavations,  blasting,  tunneling,  and  drilling. 
Such  distrubances  probably  would  not  reduce  the  spring  flow.  However, 
the  results  could  range  from  minor  to  major  redistribution  of  spring 
flow.  Major  redistribution  of  flow  could  render  the  collection  system 
ineffective. 

MANAGEMENT  OF  THE  HOT  SPRINGS 

Although  the  Hot  Springs  National  Park  is  comprised  of  many  inter- 
esting facets  such  as  panoramic  views,  hiking,  picnicing,  and  camping 
facilities,  and  a  wide  variety  of  flowers  and  trees,  the  hot  springs 
are  the  focal  point.  Proper  management  of  the  springs  will  also  preserve 
the  hydrologic  system  through  which  the  springs  work.  Measurement  of 
spring  flow,  temperature,  streamflow,  and  water  quality  will  enable  the 
Park  Service  to  monitor  any  significant  changes  in  the  system.  Proper 
management  also  depends  upon  a  knowledge  of  cause-effect  relations  be- 
tween man's  activities  and  their  effect  on  the  hot  springs,  and  authority 
to  control  man's  activities  that  can  affect  the  hot  springs. 

Spring  flow  is  now  monitored  by  a  recording  gage,  located  in  the 
reservoir  underneath  the  National  Park  Service  Headquarters,  and  by 
periodic  measurements  of  flow  of  the  two  display  springs.  Temperature 
recorders  are  operated  on  the  Fordyce  and  Quawpaw  Springs  and  manual 
temperature  measurements  of  the  display  springs  are  made  periodically. 
Two  rainfall-runoff  recording  gages  are   operated  on  the  tributaries  to 
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Hot  Springs  Creek;  one  is  located  along  Park  Avenue,  near  the  Vapors 
Club,  and  the  other  along  Whi ttington  Avenue,  near  the  Weyerhaeuser 
offices.  Nearly  all  the  individual  springs  have  been  sampled  and  the 
waters  analyzed. 

Continued  measurement  of  springs  and  streams,  collecting  and 
storing  the  data,  assembling  and  collating  the  data  annually,  and 
analysis  of  current  and  historical  data  would  permit  an  evaluation 
of  the  desirability  of  discontinuing  certain  parts  of  the  data  collec- 
tion, such  as  continuous  records  of  streamflow,  after  a  sufficient 
period  of  record  has  beeri  established.  Further,  annual  data  analysis 
would  provide  an  early  indication  of  changes  in  the  spring-flow  system 
and  possibly  could  be  used  to  devise  and  apply  remedial  measures,  if 
necessary. 

Ideally,  the  hot  springs  could  best  be  preserved  by  maintaining 
the  recharce  and  discharge  areas  of  the  springs  in  as  near  their  natural, 
presettlement  state  as  possible.  However,  multiple  use  of  the  springs 
and  the  spring  recharge  and  discharge  areas  can  be  accomplished  without 
undue  stress  on  the  spring  system.  Successful  multiple  use,  however, 
may  restrict  the  degree  and  type  of  development  in  the  spring  recharge 
and  discharge  areas. 

Developments  such  as  construction  of  impermeable  covers  over  the 
recharge  area  and  pumping  of  ground  water  could  reduce  flow  in  the  system. 
Surface-  and  subsurface-waste  disposal  in  the  recharge  area  could  cause 
contamination  or  pollution  of  the  springs.  In  the  discharge  area, 
special  precautions  may  be  required  to  seal  the  springs  against  contamina- 
tion and  to  prevent  spillage  of  liquids  or  soluble  solids  and  heavy  use 
of  pesticides  and  herbicides  on  the  hillside  above  the  springs. 
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